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In Part I of the present series of papers, the experimental data for the 
specific heat at constant volume of the four metals aluminium, copper, 
silver and lead were analysed and discussed. The method of analysis was 
a determination, with the aid of a table of Einstein’s specific heat function, 
of the frequency of atomic vibration which, if assumed to be common to 
all the atomic oscillators, would give for the specific heat of the solid its 
observed value at any given temperature. The frequency thus determined 
is termed the effective average frequency of the atomic oscillators and is 
itself a function of the temperature. When this is represented as a graph, 
it discloses a progressive diminution of the frequency in the upper part of 
the temperature range covered by the data. Considered in relation to the 
atomic vibration frequencies themselves, this diminution is least conspi- 
cuous in the case of copper, while the three other metals aluminium, silver 
and lead form a sequence in which it exhibits an increasing magnitude. 


The explanation of these results is not far to seek. The Einstein formula 
for the energy of an atomic oscillator as a function of its frequency and 
the expression for the specific heat deduced therefrom by differentiation rest 
on the assumption that the oscillators are harmonic. This assumption is 
however not strictly justifiable. The expansion of volume which all the 
metals exhibit with rise of temperature can only be understood if we recog- 
nize that the atomic oscillators which are the carriers of the thermal energy 
are anharmonic in their modes of vibration. Such anharmonicity implies 
that the vibration frequencies of a quantised oscillator do not form a 
harmonic sequence but that the successive increments of frequency fall off 
in a regular progression. Simultaneously also, the dimensions of the 
oscillator would exhibit a progressive increase which is connected with the 
magnitude of the vibration and its anharmonicity and hence also with the 
changes of frequency in a quantitative fashion. It is worthy of remark in 
this connection that the volume expansion between 0° and 300° absolute 


Al 59 


INTRODUCTION 


4 


60 


Sir C. V. RAMAN 


in the case of the four metals under consideration is least for copper, being 
1%, while the three other metals, aluminium, silver and lead form a sequence 
in which it progressively increases, being 1-2%, 1-3% and 2:2% respectively 
for these metals. Thus, there is a clear correlation between the rate of 
thermal expansion and the rate of change of the atomic vibration frequencies 
which these metals exhibit at higher temperatures as deduced from their 
specific heat data. 


As is well known, the thermal expansion coefficient of a metal tends 
to zero at the lowest temperatures and at higher temperatures increases 
pari passu with the specific heat. It follows from these facts and the 
remarks made above that the specific heats of these metals can be adequately 
represented by a summation of the Einstein functions for frequencies 
assumed to be independent of the temperature only over the limited part 
of the temperature range where the overall expansion of the solid is small 
and may be ignored without sensible error. 


In Part II of the present series of papers, the characteristic frequencies 
for the four metals were evaluated by an approximate method based on their 
relationship to their elastic constants. As was remarked in that paper, 
the specific heats thus calculated showed distinct deviations from the 
observed values in the steeply falling part of the curve where its course is 
specially sensitive to the values chosen for the frequencies. The need for 
a more precise evaluation was thereby rendered evident. In the present 
paper, it will be shown that this revision may be effected by making use of 
the specific heat data themselves in the lower part of the temperature range. 


2. THE CHARACTERISTIC FREQUENCIES 


Table I below gives the four characteristic frequencies for the four 
metals calculated approximately from their elastic constants and used in 
Part II. In the adjacent columns are shown the revised values adopted in 
the present paper. The latter have been determined by trial so as to give 
the best possible fit with the actual course of the specific heat curve in the 
lower part of the temperature range in each case. Multiplying each fre- 
quency by its degeneracy and dividing the sum total by 21, we obtain the 
average frequency of the atomic oscillators and this is also shown at the 
foot of each column. 


On an examination of the figures in Table I, it will be seen that in all 
the four cases, the frequencies as now revised exhibit certain general 
features in common. Their magnitudes follow each other in the order |, 
Ill, Il and IV. This is also the order to be expected theoretically, for the 
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TABLE I 


Characteristic frequencies in cm. 


| Copper | Silver Lead 
| 


From From , From From From From | From 
elastic i elastic | specific | elastic specific | elastic | specific 
jconstants constants} heats j|constants} heats j|constants} heats 


353 179 
177 109 
287 179 
145 70 


Average ..| 211 


i 


cubic spacings are smaller than the octahedral ones and hence should 
correspond to higher frequencies of vibration, while the longitudinal oscil- 
lations in each case should have a notably higher frequency than the 
transverse ones. Comparing the frequencies determined from the specific 
heat data with the approximate values estimated from the elastic constants, 
we also notice certain general features: the frequency of mode I as now 
deduced is slightly smaller in all the cases; the frequency of mode II is 
slightly larger in all the cases except lead; the frequency of mode IV is 
distinctly larger in all the four cases, while the frequency of mode III 
is notably diminished in all the cases except aluminium. These alterations 
have resulted in an increase by a few per cent. of the average of the frequencies 
for aluminium and silver and a diminution by a few per cent. of the average 
in the cases of copper and lead. They have also resulted in bringing the four 
frequencies for each metal into the proper relationships with each other to 
be expected in physical grounds, already remarked. 


3. THE EVALUATION OF THE SPECIFIC HEATS 


Using the values of the characteristic frequencies now adopted, the 
specific heat for each of the four metals has been calculated by summing up 
the Einstein functions for the four frequencies and adding the contribution 
from the continuous spectrum, the upper limit of frequency for the latter 
being the same as the lowest of the four characteristic frequencies. From 
the specific heats thus evaluated, the effective average frequencies of the 
atomic oscillators has also been evaluated as a function of temperature. 
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In Figs. |, 2, 3 and 4 above, these calculations have been represented 
graphically. The abscisse in the graphs are the absolute temperatures, 
while the ordinates for the specific heats and for the effective average 
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frequencies appear respectively on the left- and right-hand sides of the 
figure. It will be seen that the effective average frequency at the upper end 
of the temperature range is practically the same as the arithmetical average 
frequency shown at the foot of the columns in Table I in each case. 
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The experimentally determined values of the specific heat and the 
effective average frequencies calculated therefrom have also been shown in 
the graphs alongside of the theoretically computed curves in Figs. 1, 2, 3 
and 4. As already stated, the characteristic frequencies employed for the 
calculation were so chosen as to give the best possible fit with the specific 
heats in the lower part of the temperature range, i.e., between 15° and 150° 
absolute in each case. It will be seen that in this part of the temperature 
range, the observed values fall close to the theoretical curves. On the 
other hand, in the upper part of the temperature range, i.e., between 150° 
and 300° absolute, the specific heats observed are consistently higher than 
the calculated values. The reason for this is evident from the plots of the 
effective average frequency of atomic vibration. In all the four cases, the 
atomic vibration frequencies actually effective in the upper part of the 
temperature range are no longer those deduced from the specific heats 
at low temperatures but show a progressive diminution with rise of tempe- 
rature, for the reasons already explained in the Introduction. 


SUMMARY 


The specific heat data themselves enable us to evaluate more precisely 
the four characteristic frequencies of each metal determined approximately 
in Part If from their elastic constants. Making use of the new values, the 
complete specific heat curve is theoretically deduced and shows good agree- 
ment with the observations in the lower part of the temperature range, but 
deviates observably in its upper part, as is to be expected in view of the 
progressive fall in the frequency of the atomic oscillators with rise of tem- 
perature indicated by the analysis of the data in Part I. This diminution of 
the frequencies is a consequence of the anharmonicity of the oscillators 
which also results in the thermal expansion of the metal. Copper, aluminium, 
silver and lead form a sequence in the order of increasing coefficients of 
thermal expansion as also in the magnitude of the temperature coefficients 
of atomic vibration frequency indicated by their specific heat data. 
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OXIDATIONS WITH SILVER IODIDE DIBENZOATE 


Part II. Preparation of s-Azo Compounds from Some Aromatic Primary Amines 


By P. S. RAMAN 
(Department of Organic Chemistry, University of Madras, Madras-25) 
Received December 18, 1956 
(Communicated by Dr. K. N. Menon, F.A.sc.) 


IN a recent paper! from this laboratory the use of the Sinomini complex 
formed from iodine and silver benzoate for the scission of a-glycols was 
reported. The present communication deals with the oxidation of some 
substituted anilines. Aromatic primary amines are rapidly oxidised by the 
reagent with the production of deeply coloured products. A detailed investi- 
gation of the action of this reagent on p-toluidine, p-chloraniline, p-brom- 
aniline, p-nitroaniline and 2-chloro-4-nitro-aniline has shown that the main 
product of the oxidation in each case is the corresponding symmetrical 
azo compound, when molar proportions of silver iodide dibenzoate and 
the aromatic primary amine are employed. The reaction thus proceeds 
according to the general equation 


2R-NH, -++ 2 (C,H;-COO), Agi ——> R-N = N-R + 2Ag 1 
+ 4 C,H,-COOH. 


Silver iodo dibenzoate thus constitutes another useful member of a 
group of oxidising agents like permanganate,” ferricyanide,* chromic acid 
anhydride, sodium hypobromite,® phenyl iodosoacetate® and lead tetra- 
acetate?» § for the oxidation of aromatic primary amines to the corresponding 
s-azo compounds. 


The reaction probably involves a free radical mechanism and proceeds 
through the highly reactive intermediate radical R-NH° which undergoes 
dimerisation to the hydrazo compound followed by further oxidation to 
the azo compound. 


(CgHs—COO)2 Agi 
2 R-NH, — 2 R-NH* ———» R-NH-NH-R 


R-N = N-R. 


The intermediate R-NH° may enter into side reactions and _ this 
probably accounts for the formation of dark coloured products encounter- 
ed in the oxidations. This mechanism finds support in the work of Pausacker 
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and others on the oxidation of aromatic primary amines to azo compounds 
with pheny! iodosoacetate® and lead tetra-acetate? for which analogous re- 
action mechanisms have been proposed by them. The mechanism is also 
partially supported by the fact that hydrazo benzene is very rapidly 
oxidised to azo benzene in very high yields (94%) by the action of silver 
iodide dibenzoate. But aniline itself, on treatment with the complex under 
mild conditions, gave a dark, gummy product from which only traces of 
azo benzene could be isolated. 


Although the conditions for the maximum yield of the azo compounds 
were not investigated, the yields obtained in the different cases were fairly 
good, excepting in the case of p-toluidine. The oxidation of p-toluidine 
was very rapid whereas it was relatively slower in the case of p-nitroaniline 
and still slower in the case of 2-chloro-4-nitroaniline. The yields of azo 
compounds were somewaht widely different in the case of the different amines; 
those amines containing electron-attracting groups tend to give relatively 
higher yields of azo compounds. Thus p-toluidine gave 28%, p-nitroaniline 
44% and 2-chloro-4-nitroaniline 56% of the corresponding azo compound. 
Benzoic acid, arising from the oxidising agent, was obtained in high yields 
in all cases after the oxidation. 


EXPERIMENTAL 


Silver iodide dibenzoate was prepared by warming 9-2 g. (0-04 mole) 
of silver benzoate and 5 g. (0:04 atom) of iodine in 140 c.c. of dry benzene, 
with exclusion of moisture. This was cooled to room temperature and the 
amine (0:02 mole) added either as such or in benzene solution. The mix- 
ture in almost all the cases became warm and was afterwards shaken at room 
temperature and then heated on the water-bath for a short time. By this 
time the bulky silver iodide dibenzoate got replaced by the heavier silver 
iodide and the benzene solution became deeply coloured. The hot benzene 
solution was filtered, the residual silver iodide washed with warm benzene 
and the combined benzene filtrate worked up for the isolation of the azo 
derivatives as indicated under the different compounds. The identity of 
the difierent azo compounds was established by mixed m.p. determinations 
with authentic specimens prepared according to the general method of 
Bogoslovskii.® 


4, 4’-Dimethyl azo benzene.—After adding the p-toluidine dissolved in 
dry benzene to the complex, the mixture was shaken for 10 mins. and then 
heated on the water-bath for about 15 mins. The colour of the benzene 
solution which was brownish yellow at first turned dark brown by this time. 
The benzene solution from the reaction was concentrated to small bulk 
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under reduced pressure, cooled and the crop consisting mostly of benzoic 
acid was filtered off. The benzene filtrate was brought to dryness and the 
residue triturated with alcohol. The portion insoluble in alcohol was 
brown in colour and melted at 135-40°. Wt. 0-6 g. (28%). Crystallisation 
from alcohol gave pure 4, 4’-dimethyl azo benzene, yellow crystals, m.p. 
and mixed m.p. 143-44°. 


4, 4'-Dichloro azo benzene.—p-Choraniline in benzene solution was added 
to the silver iodide dibenzoate complex, mixture shaken for 15 mins. and 
then heated on the water-bath for 10 mins. Concentration of the dark 
brown benzene filtrate under reduced pressure to about 40c.c. and cooling 
in ice gave the azo compound mixed with benzoic acid. The azo com- 
pound (yield 1 g.; 40%; m.p. 180-85°) separated by bicarbonate, crystal- 
lised from alcohol and then acetone (charcoal) gave yellow needles, m.p. 
and mixed m.p. 184-85°. 


4, 4'-Dibromo azo benzene.—The oxidation of p-bromaniline was done 
exactly as in the case of p-chloraniline. The colour of the benzene solution 
which was brown initially changed to almost dark towards the end. The 
benzene filtrate, worked up as above, gave | -4 g. (41%) of the azo compound 
as a yellowish brown solid, m.p. 197-200°. Crystallisation from acetone 
(charcoal) gave yellow needles, m.p. and mixed m.p. 203-05°. 


4, 4’-Dinitro azo benzene.—Powdered p-nitroaniline was added to the 
complex and the mixture heated on the water-bath for 75 mins. The deep 
red benzene filtrate on concentration to small bulk, gave the azo compound 
mixed with benzoic acid. Separation by bicarbonate gave the azo compound 
(1-2 g.; 44%) which on crystallisation from acetone (charcoal) gave orange- 
red needles, m.p. 223-24°. 


2, 2'-Dichloro-4, 4'-dinitro azo benzene.—After addition of 2-chloro- 
4-nitroaniline to the complex, the mixture was very gently boiled on the water- 
bath for 14 hrs. The benzene solution was brought to dryness under reduced 
pressure, residue treated with bicarbonate solution and then triturated with 
cold alcohol. Filtration gave a brown product (1:9g.; 56%) melting at 
about 270°. Crystallisation from benzene gave the azo compound as fine, 
brownish red crystals, m.p. and mixed m.p. 275-76°. 


Azo benzene from hydrazo benzene.—Hydrazo benzene (3-7 g.) was added 
to the complex prepared as above and after shaking the warm mixture for 
15 mins., it was heated on the water-bath for 10 mins. The colour of the 
benzene layer which immediately changed to yellow, turned reddish yellow 
and finally reddish brown. The filtered benzene solution was brought to 
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dryness, the residue taken up in ether, washed with bicarbonate solution 
and then water, the ether solution dried with anhydrous sodium sulphate 
and the solvent removed. The residue on a single crystallisation from dil. 


alcohol gave 3-45 g. of azo benzene as orange-red crystals, m.p. and mixed 
m.p. 67-68°. 
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SUMMARY 


The formation of s-azo compounds by the oxidation of some aromatic 
primary amines by silver iodide dibenzoate is described. The reaction 
probably proceeds through a free radical mechanism in which the highly 
reactive R-NH?® radical plays a part. 
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By S. K. SRINIVASAN AND N. R. RANGANATHAN 
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Received December 3, 1956 
(Communicated by Dr. Alladi Ramakrishnan, F.A.sc.) 


ABSTRACT 


In the new approach to the cascade theory suggested by Alladi Rama- 
krishnan and one of us (S. K. S.), we deal with the number of electrons 
produced between 0 and ¢ in a soft cascade, each of these electrons having 
an energy greater than E at the point of its production as contrasted 
with the usual approach where we are interested in the number of electrons 
with energy greater than E at ¢. We present here numerical calculations 
of the mean numbers on the basis of the new approach. 


IN a previous paper (1956) in these proceedings (hereafter referred to as 
Paper 1), Ramakrishnan and one of us (S. K. S.) suggested a new approach 
to cascade theory which yields equations having elegant asymptotic 
solutions. In the normal approach we are interested in n(E; t) the number 
of particles above a certain energy E at a particular thickness ¢ in a shower 
initiated by a particle of known energy or energy distribution. On the other 
hand, we now ask for N(E; ¢) the number of particles produced between 
0 and ¢ each with energy greater than E at the point of its production. It 
is well known that n(E; 1) and N (E; 2) are stochastic variates and a com- 
prehensive treatment should deal with their probability distribution func- 
tions. A more limited and tractable problem is to obtain the first and second 
moments of the distribution. In the case of electron-photon showers the 
first and second moments of n (E; t) have been tabulated by other workers.* 
We now here present the numerical calculations relating to the first moment 
of N(E; 


In Paper I, Ramakrishnan and one of us (S. K. S.) have proved that 
the mean number of electrons produced between 0 and ¢ with their primitive 
energies above E, in a shower initiated by a single electron of energy E, is 
given by 


* For a comprehensive account of the treatment, see Ramakrishnan and Mathews (1954). 
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e{N(E|E,; 9} 


where 
l 4 l 
Bs = 2 +1) (s+ »}- 
4 
Cs = s+] + G «) 5G — (3) 
A; =4(As + D) — ${(As — D)? + 4 B,C}. (4) 
Hts = 4(As + D) + 4 {(As — D)* + 4 (5) 
As =(5 + «){ log +y + st +4 (6) 
D =) a= 0246. (7) 


For fairly large thicknesses (say ¢ > 4), we can neglect the term contain- 
ing e~“s* and (1) can be written as 


g+ico 
B,C, E,\*" 
¢+ico 
As As(s — I) E) 
ico 
= fen ds ds (8) 
@¢—ico o@-—ico 
where 
B,C E 
$1 (S) = log — log(s — 1) + y(s — D3» = log (9) 
$2(S; t) = log Goto — log(s — 1) + r(s— 1) — tA, (10) 
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Noting that ¢{N (E/E; f)} is a function only of E,/E or y = log E,/E we 
can write it as «e{N(y;2)}. We next observe that as t—» oo the second in- 
tergral tends to 0 and hence we have 


¢+ico 


T—ico 


As s-—>2, As—0O and B,C,/u, remains finite and non zero at s = 2. 


Hence ¢,(s) and ¢.(s; f)-—»-+ coas s—»2. From (2) to (7), it follows 
that for very large real part of s, 


As ~(F+ a) (log s + +4 


Bs =~ 
B;C 2 
+4) s(s + 1) logs 
B,C 2 
+2) s(s + 1) logs 
This 
and 


tend to — oo logarithmically. Since y(s — 1} -»-+ oo much more rapidly, 
¢, (s) and ¢.(s; ft) tend to + co, as s—»co. Hence ¢, and ¢, must have 
a minimum as s increases along the real axis from 2 to oo. Let s, and s, 
be the points at which 


B, C 
log — log’'s — 1) — ta, 
(44) =0, (2) (13) 
ds J sus, OS < 
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We can shift the line of integration in each of the integrals so that it 
passes through the saddle points. This can always be achieved since the 
only restriction on the line of integration is that it should be to the right of 
imaginary axis and also to the right of all singularities, i.e..o>2. Having 
chosen the contour to pass through the saddle point, we have the well-known 
approximation of the saddle point method obtained by replacing each of the 


functions ¢, and ¢, by the first three terms of its Taylor expansion about the 
saddle point. 


Thus 


= — (14) 


(83 1) = — (8) (15) 


SeS, 
where 


i€é = s— s, in (14) and in = s — Sy in (15). 
Substituting (14) and (15) in (8), we have 


+00 


+00 
= fe {8,5 dy (16) 


7 
(primes indicate differentiation with respect to s). 
Hence 
N(y3} = e%(5,) t) (17) 

The mean numbers « {N (y; f)} for various values of y and t have been 
calculated numerically using (17). In order to obtain the saddle points, 
it was found necessary to sub-tabulate at intervals of -025, the basic func- 
tions As, Bs, Cs, As and ws tabulated by Janossy and Messel (1951) at in- 
tervals of -1. The mean number of electrons produced in infinite thickness 
and fairly large thickness are given in the following table. For convenience 
of comparison, we have also tabulated « {n(y;1)}, the mean number of 
electrons that exist at ¢ with an energy above E.f 


+ These are taken from the table of Ramakrishnan and Mathews (1954), 


— co 
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e {N (y; 
«{n(y; t)} is given in brackets. 


yit 4 6 8 16 oo 
8 293 528 851 1372 1406 
(94-2) (149) (146) 
7 124 249 355 502 508 
(48-8) (61-3) (49-6) 
6 64:5 111 146 183 184-5 
(23-4) (22-9) (15-1) 
5 30-6 41-6 57-9 66-4 66-5 
(10-2) (7-55) (4-01) 
4 10-9 14-7 18-9 20:51 20:53 
(3-85) (2-11) (1-897) 
3 6-65 7-86 8-56 8-82 8-82 
(1-21) -475) (- 160) 


(y; unlike {n(y; has an interesting asymptotic property. 
While « {n(y; t)}} +0 as t+ 00, « {N(y; 4} tends to a finite value for any 
finite y as f—» oo. The smaller the y we observe from the table, the more 
quickly does « {N (y; ¢)} approach the limit. For almost all the y given in 
the table the limit is nearly reached at ¢ = 16. 


For small thickness, we cannot neglect the term containing e-* in (1). 
Then the saddle point formula as adapted here cannot be applied to this 
term. An alternative method of calculation of (1) for small thickness is being 
attempted and the mean numbers for small thickness will be reported in 
a subsequent contribution. 


We are deeply indebted to Dr. Alladi Ramakrishnan for suggesting this 
problem. 
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ON THE DETERMINATION OF THE SLOWING 
DOWN LENGTH IN ORDINARY WATER USING 
A PULSED NEUTRON SOURCE 


By R. RAMANNA AND M. P. NAVALKAR 
(Atomic Energy Establishment, Bombay) 


Received November 22, 1956 
(Communicated by Prof. B. Peters, F.A.sc.) 


INTRODUCTION 


A METHOD of determining the slowing down age in ordinary water from a 
fast neutron source to thermal energy, using a pulsed neutron source is de- 
scribed. The fast neutrons were produced from the D-Be reaction using a 
cascade generator. This reaction emits neutrons of maximum energy of 
4-5 Mev with a peak at about 3-5 Mev. Pulsed neutron methods have been 
used by earlier workers in the measurement of the mean life of thermal 
neutrons in water (Haworth ef al., 1942),1 (Von Dardel et al., 1953),* 
(Antonov ef al., 1955),* the diffusion length of ordinary water and other 
substances (Von Dardel, 1953),* the diffusion parameters of ordinary water 
(Von Dardel er al., 1954)° and the diffusion and the slowing down constants 
of ordinary water and beryllium oxide (Ramanna er al., 1956).® 


THEORY 


The slowing down of neutrons in a moderating medium can be ex- 
plained either on the Age Theory or Multi-Group Theory. The Age Theory 
assumes ‘a continuous process of energy loss of source neutrons to thermal 
energy whereas in the Multi-Group Theory the energy interval from source 
to thermal is divided into finite number and there is no energy degradation 
of neutrons in any group until they go over to the next group. The Multi- 
Group Theory approximates to the Age Theory when the number of energy 
intervals tend to infinity. But in the case of light moderators like water 
where the number of collisions to thermalise is small, the Age Theory is a 
poor approximation to the slowing down process. 


The spatial and time distribution of thermal neutrons from a pulsed 


source of fast neutrons has been given on the two theories as follows (Mani 
and Iyengar, 1957)’ 


Age Theory 


S 
ps (rt) = 
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exp [- | 
x 4 {7 + (t Tef)/Tes} J (1) 
[7>) + Ls? (t — 


where 
tcf = Slowing down time of source neutrons to thermal energy. 
Tes = Mean life of thermal neutrons in the moderator. 
7) = Age of source neutrons to thermal energy. 
L, = Diffusion length of thermal neutrons in the moderator. 
S) = Source strength. 
Two Group Theory 


—t 
TesSo 


ps (r,t) = [er I (as, Vu) — eter I (as, vw) 


= erf(x) = 


where the subscript ‘ s’ denotes the slow term and ‘ f’ denotes the fast term. 


a 
I(a, \/n) we” e¥” upto a time 1, after which it decreases rapidly. Since 
Tef Kes after a time te, I (af, Wu) <I (as, and can be neglected. 
Hence between and being less than 
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and so 
[1 — — — (5) 
Unlike the Two Group Theory, the Age Theory expression does not show 
such an exponential character in any interval of time. 
Also 
d rps(r,t) 
6) 
But 
— 


Tes — Tef 


p= (7) 


(since te¢ < Tes) 


where Ly* = Square of the slowing down length of source neutrons to 
thermal energy. 


Therefore, L¢ can be determined using the above expression from a measure- 
ment of the spatial distribution of neutron from a pulsed source. 


EXPERIMENTAL ARRANGEMENT 


The neutron pulsing system.—The fast neutron bursts were produced by 
pulsing the extraction electrode of the radio frequency ion source of the 
cascade generator. The pulse width was kept at 19 y-sec. and the 
repetition rate at 560 cycles per second. 


Detector.—The thermal neutrons were detected by means of a lithium 
iodide scintillation counter. The crystal (1-25 cm. in diameter, | cm. long) 
was mounted at the end of a long perspex rod (3 ft. long, | in. in diameter) 
which was attached to the photo cathode of a 6262 E.M.I. photo multiplier. 
The pulse from the disintegration of lithium due to thermal neutrons. is 
much higher than the background radiation and can thus be easily discrimi- 
nated. The light guide assembly was provided with a spigot at its end so 
that the position of the crystal could be fixed at known distances and in the 
same plane from the centre of the target. This was done by inserting the spigot 
into various positions on a perspex strip attached to the target holder. The 
strip was provided with holes at fixed positions to take in this spigot and was 
so arranged that it was kept at the centre of the assembly. A rectangular 
tank (60 cms. cube) of perspex was filled with approx. 55 gallons of dis- 
tilled water. The sides and bottom of the tank were covered with -5 mm. 


thicl 
the 


{in ona 
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thick cadmium on the outside to absorb thermal neutrons entering or leaving 
the tank. The block diagram of the experimental set-up is given in Fig. 1. 
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Block Diagram of the Experimental Apparatus 


The Time Analyser—The thermal neutrons were analysed in time by 
a 10 channel time analyser (lyengar, 1956). For the experiment each channel 
width was adjusted to 54 w-sec. with an initial delay of 200 p-sec. 
In order to get more points in between, the initial decay was shifted to 227 
u-sec. keeping the channel width same. The two sets of readings were 
normalised by the thermal neutron counts taken in the same interval of 
time for both the delays. In this case the interval of time was 0-150 
psec. 


Experimental neutron distribution for some of the detector positions 
is shown in Fig. 2a and 26 while the theoretical curves on both the Age 
and the Two Group Theory for r = 30cms. are given in Fig. 2c. It is 
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F 22.5 

sk A = -003826+ .000017 secs” 
~ 2F 
as 

° 

TIME IN SECS. 
Fic. 2a 


noticed that in the region between 200 and SOO yp-sec. the experi- 
mental distribution plotted on a log scale is a straight line when the detector 
is at 22-5 cms. from the source. The straight line portion increases in time 
duration as the detector is moved away from the source. The values of the 
decay constants at various detector positions from the source have been 
determined by least square fitting and are given in Table I. 


Calculations of the decay constant were made on both Age and Two 
Group Theories. The constants used in the calculations were as follows: 
Using Marshak’s formula for the age of source neutrons to indium energy 
and correcting it for thermal energy, the value obtained for age is about 
39 cms.” for 3 Mev. and 44cm.? for 3-5 Mev. mono-energetic neutrons. 
Taking a value of about 41 cms. for age, 7¢; = 216 u-seconds,® tz = 10p- 
seconds,® L, = 2:88 cms.® and using expression (5), the Two Group Theory 
gives a decay constant of 0-0037 yu-sec.—! which is in better agreement with the 
experimental value. Unlike the Two Group Theory, the Age Theory gives 
different decay constants at different detector positions while those obtained 
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TABLE | 


Distance r from Time in which A 


the source is calculated Decay constant A Error in A 
in cms. in p-sec. p-sec.-! dx 
20 227-389 0-003558 + 0-000015 
227-497 0-003526 -- 0-000017 
25 227-497 0-003450 -- 0-000019 
227-497 0-003455 0:000024 
30 335-686 0-003434 0-000022 
335-686 0-003429 -+- 0-000024 
Vos 30Cms. 
0:003434 + 0000028 secs” 
EXPERIMENTAL CURVE 
ce 
th 
TIME IN J SECS. 
Fic. 2 b 


experimentally are independent of the detector positions as long as they 
are greater than 22-5 cms., thus showing that the Two Group Theory 
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THEORETICAL CURVE ON AGE THEORY 
FOR ¥ = 30 cms. 


THEOCRETICAL CURVE ON TWO GROUP THEORY 
FOR ¥* 30CmMms. 


LOG. p 


2 4 


° 250 500 75° 1000 


TIME IN M Secs 


Fic. 2c 


applies better than Age Theory in case of water. The disagreement of Age 
Theory is to be expected since it is well known that the Age Theory is in- 
applicable in case of light moderators like water. The experimental value 


of the decay constant is in agreement with the value obtained previously, 
(Ramanna et al., 1956).® 


The values of the decay constants at detector distances of 10 and 15 cms. 
from the source are much higher than expected on the Two Group Theory. 
This may be explained on the basis that at such close distances there is a 
large concentration of fast neutrons which decay more rapidly than thermal 
neutrons and hence give higher value of the decay constant. The decay 
constant therefore decreases with distance from the source and _ then 
remains constant. 


Fig. 3 shows log — “g-—- against * and is seen to be a straight line. 
0 


The slope of this line at various times is given in Table II. Using the 
average value of these slopes and taking tes = 216 p-sec,® 7-7 = 10 micro-sec.® 


- 
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S\.OWING OCOWN CENGTH FOR SOURCE 
NEUTRONS TO THERMAL ENERGY 


SLOPE = .15334 + .00022 cm' 


« 


S(v, t) 


oO 25 5 75 10 12-5 168 17-5 20 22-5 25 27-5 30 32-5 35 


IN Cms. 
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TABLE Il 
Time in preg ty 
microseconds log s "’—r Error in the slope 
362 0- 15307 t+ 0-00021 
416 0- 15334 + 0-00022 
470 0- 15297 + 0-00022 


and L, = 2-88 cms. in expression (7), we get the age of source neutrons 
to the thermal energy as 41-05 cms.” The value of the age changes to 
39-47 cms.” when vez is taken as 20 p-sec. If we use the approximation in 
expression (7), i.e., te = 0, we obtain 42-65 cms.* as the value for the age. 
It can be noticed that in altering the value of zs, the value of the age does 
not change more than + 4-0%. The statistical error due to counting is 
very small compared to the above error and can be neglected. It is 
unlikely that the value quoted above for age 41-05 cms.? has been deter- 
mined to better than + 4:0%, taking into account both the approxi- 
mations involved in theory and the experimental errors introduced due 
to the finite dimensions of the detector, the source and the moderator 
assembly. It can, therefore, be concluded that the experimental value 
obtained for age by this method is consistent with the value expected for 
the D-Be neutrons. 


We wish to express our thanks to Mr. S. B. D. Iyengar, Mr. G. S. Mani 
and Mr. P. K. Iyengar for helpful discussions. Our thanks are also due 
to Mr. S. K. Ambardekar for running the accelerator and Mr. F. R. Bhathena 
for his help in taking readings with the time analyser. 
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1. INTRODUCTION 


Let A = (ays) be an arbitrary square matrix of order n and c (A) acharacter- 
istic root of A. Hirsch® proved that 


|c(A)| <n. max. | ays |, (1+1) 


and Browne® obtained that - 
|c (A) | <(R+ T)/2, (1 +2) 
where 


n n 
R = max. R, = max. 3 | ars |, T = max. T, = max. J | dps j. 
(r) (r) OF1 (s) (s) 


Further, Browne? showed that 
s<|ec(A)|? <G, (1-3) 
where s and G are the smallest and greatest characteristic roots of AA’. 


(1-1), (1-2) and other results established by several authors (for a list 
of references see [1]) give upper limits to | c (A) | in terms of the elements of 
A, while (1-3) gives limits to | c(A)| in terms of the characteristic roots 
of AA’. In the present note an upper limit for | c(A)| has been found in 
terms of the characteristic roots of the associated Hermitian matrices 
(A + A’)/2 and (A — A’)/2i, and this result has been generalised for the 
characteristic roots of the product matrix AB of two n-square real or complex 
matrices A and B. An upper limit for | c(AB){| has also been found in 
terms of the elements of A and B. This result contains (1-2) as a particular 
case. 


In establishing these results we use the well-known facts, see [4; p. 149], 
concerning the operator norm 


|| A || = max. || Ay || for | y || = 1, when || y | = Vy. 
48 
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In what follows R(A) = max. R, (A) = max. | ars | and T(A) = 


s=1 


max. Ts (A) = max. x | dps |. 
(8) (8s) 


2. BOUND THEOREM FOR c(A) 


THEOREM 1.—Let c(A) be an arbitrary characteristic root of A. Then 
A A’ — A’ 
) + max. | le (2.1) 


Proof.—Any square matrix A = (A + A’)/2 + (A — A’)/2 = M+ iN, 
say, where M and N are Hermitian. Corresponding to an arbitrary 
characteristic root A of A, there exists a non-zero column vector x, such that 
Ax = Ax. Then 


| c (A) | < max. |c 


[A]. = = Axl 
= ||(M + iN) x|| 
<|M+iN||. 


where for any square matrix B, || B ||? = max. || By || = max. (By, By) 


= max. (j'B’By) = = (BB’), denoting the greatest 7 


characteristic root, and in particular if B is Hermitian || B || = max. | c (B) |. 


Since 
= > 0, 
we have 
|A| <||M-+iN || 
<||M||+1N]] 


= max. | c (* +) + max. 


A- )|- 
This establishes (2-1). 


It may be observed that (1-3) and (2-1) give the same upper limit for 
|c(A)| when A is Hermitian or skew-Hermitian. There are some other 


matrices also, for example (' 3; Re 1) , which are not of some special type 


but for which (1-3) and (2-1) give the same upper limit for | c (A) |. 
limit for | c (A) | given by (2-1) is not better than that given by (1-3), but, 


| 
=) 
er 
2. 
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‘ 


86 N. A. KHAN 
in general, it is easier to find the characteristic roots of (A + A’)/2 and 
(A — A’)/2i than those of AA’. 

3. BOUND THEOREMS FOR c (AB) 


We now consider two theorems which give upper limits to the absolute 
value of an arbitrary characteristic root of the product matrix AB. 


THEOREM 2.—Let c(AB) denote an arbitrary characteristic root of AB. 
Then 


| c(AB) | < | max. + max.|c 


x (max. + max c (3.1) 


Proof.—Any square matrix A = (A + A’)/2 + (A — A’)/2=M-+4+ iN, 


say, and B = (B + B)/2 + (B — B)/2 = R + iS, say, where M, N, R and 
S are Hermitian matrices. Thus AB =(M + iN)(R + iS). 


If » is an arbitrary characteristic root of AB, there exists a non-zero 
column vector x, such that wx = ABx so that 


= || ABx || 
< || (M + iN) ||. || (R + iS) || . ix 
Since 
\| x || > 0, 
we have 


= (max. | c (M) | + max. | c (N) |) (max. | ¢ (R)| 
+ max. | c (S))}) 
since M, N, R and S are Hermitian matrices. 


This establishes (3.1). It is easy to see that (3.1) can be generalised 
to the case of the product of any finite number of matrices” A,, As, 
Ax of the same order. The result in the generalised form is 


A, + A,’ 


<I max. c 


(3.2) 
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THEOREM 3.—Let c(AB) denote an arbitrary characteristic root of AB. 


Then 


Proof.—If » is an arbitrary characteristic root of AB, there exists a 
non-zero column vector x, such that px = ABx. 


This gives, as we have seen in the proof of the previous theorem, 
< R(A) + T(A),, R(B) + TO) 
>’ 


tor, by a known result [3; p. 708] 


This completes the proof of Theorem 3. It may be observed that 
(3.3) can be generalised to the case of the product of a finite number of 
matrices A,, Ag,.....-A, of the same order. The result in the generalised 
form is 


4. PARTICULAR CASES OF THEOREM 3 
(i) If we put B=TI, (3.3) reduces to 
A 
A), 


a result due to Browne,*? quoted earlier. 


(ii) If A and B are Hermitian or skew-Hermitian matrices, then R, (A) 
= T,(A) and R,(B) = T; (B) for r= 1, 2,...., which means R (A) 
=T(A) and R(B)=T(B). Hence for Hermitian or skew-Hermitian 
matrices A and B 


| c (AB) | < R(A) R@B). (3.5) 


In conclusion, | am thankful to Professor S. M. Shah for his help in the 
preparation of this paper. 
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IMPROVEMENTS IN COLORIMETRIC 
DETERMINATIONS 


Ill. Comparison of Colorimetric and Electrometric Methods for the 
Determination of Soil Reaction 


By K. L. KHANNA, F.A.Sc., P. B. BHATTACHARYA AND K. PRASAD 


(Central Sugarcane Research Station, Pusa, Bihar) 


Received October 29, 1956 
I. INTRODUCTION 


MASON AND OBENSTIAIN (1938) have, with proper manipulation, found 
that results obtained by Kuhn’s (1930) colorimetric method of pH measure- 
ment agree closely with those of the electrometric method (1931). Errors 
in measurement may, however, occur through contamination, use of 
glassware of inferior quality and a high concentration of soluble salts in 
the soil. While carrying out pH determinations of a large number of samples 
collected from Narkatiaganj area (representing the non-calcareous soil 
division of North Bihar, with a low content of soluble salts), the authors 
had an opportunity of testing the colorimetric method against the standard 
electrometric method. The results of these studies have been reported in 
the present paper. 


II. EXPERIMENTAL 


The electrometric method.—10 grams of air-dry soil are taken with 25 c.c. 
of CO,-free distilled water, stirred and examined in the Cambridge pH 
meter in the usual way. 


Colorimetric method.—A layer of barium sulphate (Pro. Rontgen) 
1 cm. thick is placed in the bottom of a dry test tube. A layer of soil about 
3cm. deep is then added and the tube filled with CO,-free distilled water to 
a depth of 9cm. Suitable indicator solutions depending upon the expected 
pH range are added and the test tube filled with water to a total depth of 
15cm. It is then closed with a paraffined cork, shaken vigorously for half 
a minute and left to settle. The colour deveioped is compared against 
Clarke’s (1928) colour chart and the reading noted. The indicators used 
in the present case were Bromo Thymol Blue and Cresol Red prepared 
according to Clarke’s procedure (loc. cit). 


Discussion 


The electrometric and colorimetric pH values of 300 samples were found 
to be in close agreement within a limit of 0-2 units in most cases. On 
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statistical examination of the data, the intra-class correlation coefficient worked 
out to 0-9879, indicating a very close agreement between the colorimetric 
and electrometric pH values. It may therefore be concluded that the 
colorimetric method may be employed satisfactorily in the estimation of 
pH in the field in case of soils with low soluble salt contents (below 150 mg,/ 
100 gm. soil) such as the non-calcareous soils of North Bihar (Prasad et al., 
1951). 
TV. SUMMARY 


pH values of 300 soil samples from Narkatiaganj area in North Bihar 
(non-calcareous soils with low soluble salt content) have been determined 
by the electrometric and colorimetric methods. The results are in close 
agreement, thus showing that the colorimetric method may be used for pH 
determinations in the field in case of soils of low soluble salt content. 
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INTRODUCTION 


THE composition of nickel ferrocyanide complex has been the subject of many 
investigations. Werner!.? studied for the first time the composition of nickel 
ferrocyanide complex and showed that Ni** got precipitated with potassium 
ferrocyanide giving a stable complex of the composition Ni, [Fe (CN)g]. 
Gaspar? and co-workers considered that the nickel ferrocyanide 
was a double complex or associated salt. Britton and Dodd‘ observed that 
nickel reacted with K,Fe (CN), to form normal salt which was converted 
by further precipitation into Ni, [Fe (CN),]. x K, [Fe (CN),], where x was 
less than one. Thermometric siudy® of precipitation of nickel ferrocyanide 
was shown to have a composition either NiK, [Fe (CN),] or Ni,K4 [Fe (CN)g];, 
the composition varying with the conditions of precipitation. X-ray dif- 
fraction studies® of nickel ferrocyanide revealed that as much as 0-3-0-4 
mole K,Fe (CN), existed with one molecule of Ni, [Fe(CN),]. It was 
assumed that the ferrocyanide was adsorbed on the surface of highly dis- 
persed crystals of Ni, [Fe (CN),] and not combined to form definite double 
salts. The present communication reports polarographic investigation of 
the nickel ferrocyanide complex. This method appeared desirable: Potas- 
sium ferrocyanide is not reduced’ at the d.m.e. and the method, therefore, 
js void of complications. Further, in a similar controversial circumstance, 
viz., investigation of the composition of the cadmium ferrocyanide 
polarographic method proved to be _ successful.® 


EXPERIMENTAL 


Nickel sulphate used was B.D.H. Analar sample. This was recrystal- 
lised in twice distilled water and dried at 50°C. for 48 hours. The exact 
concentration of the nickel solution prepared was estimated gravimetrically 
using dimethyl glyoxime. 

Potassium ferrocyanide used was Merck sample. KCl employed was 
of B.D.H. Analar quality. 


The polarograph used was assembled in our Laboratory; the electrical 
Circuit employed was similar to the one described in an earlier communica- 
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tion.® The diffusion current was obtained by measuring the potential dif- 
ference developed across a known resistance introduced serially in the cir- 
cuit; the accuracy of the measurement was adjusted to 0-1 uw amp. 


The capillary employed had the following characteristics: 
m = 0-665 mgm./sec. 
t = 4-14 sec. 
These values were recorded in 0-1 M KCI solution in open circuit, 
The method employed for investigating the nickel ferrocyanide 
complex involved the amperometric titration of a known volume of Ni+ 
solution of a fixed strength against standard potassium ferrocyanide. 
1M KCI solution was used as supporting electrolyte and 0-02 per cent. 
solution of gelatin was used as the maximum suppressor. The titrations were 
carried out at — 1-9 volts where the diffusion current attained a limiting value; 
and in an inert atmosphere of nitrogen to avoid the possible oxidation of 
nickel ferrocyanide complex,!® leading to a change in colour of the precipitate 


from green to yellowish brown; it may be added that in the present studies 
no such change in the colour has been noticed. 


RESULTS AND DISCUSSION 


Figure | gives a typical series of results on the variation of diffusion 
current during the progress of the amperometric titration of Ni‘* with 
K,Fe (CN), in presence of | M KCI used as supporting electrolyte. The 
data in Fig. 1 show that sharp end points could be obtained during the 
70] 
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Fic. 1. Amperometric titration of Ni against potassium ferrocynide. 


Concentration of ferrocynide 100 mM. Curves | to 8 refer respectively to 5, 7-5, 10-0, 15:0, 
20-0, 25-0, 30-0, 40-0 millimolar of Ni*+, applied voltage — 1-9 volts. 
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titration of Ni*+ with ferrocyanide polarographically. The observed 
quantities of potassium ferrocyanide for different concentrations of Ni+* 
are given in Table I. It is interesting to note from the data in columns 3 
and 4 (Table I) that the composition of nickel ferrocyanide complex 
corresponded to Ni;K, [Fe (CNg]4;_ this appears to follow the stoicheo- 
metric relationship according to the equation 


5 NiSO, + 4 K,Fe (CN), Ni;Kg [Fe (CN)e]4 5 K,SO,. 


TABLE 


Polarographic investigation of nickel ferrocyanide complex 


Volume of Ni*+ solution in the cell 25-00 ml. 


Volume required 
Concentration Concentration Observed on the basis of 
of Ni** solution of K,Fe (CN), titre value the formula 
in mM in mM in ml, Ni;K,[Fe(CN)g], 
(1) (2) (3) (4) 


nN 
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The results showed conclusively that over a wide range of concentration of 
Ni+ and K,Fe (CN),, the composition of the complex, as derived from the 
end points of amperometric titration curves, appeared to be Ni;K, [Fe (CN)g], 
or 3 NiK, [Fe (CN),], Ni, [Fe (CN)g]. 
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The data further show that the method could be employed with ease 
for quantitative estimation of Nit+. The experimental error involved is 
less than 1 percent. 


A detailed examination of the titration curves shows a significant 
feature in the variation of the diffusion current at the initial stages of the 
titration. While the diffusion current was expected to decrease continuously 
with the addition of the titrant leading to the familiar L-shaped curves, in 
the present investigation the current values did not decrease in the initial 
stages and remained almost unaffected; this feature was more apparent 
at higher concentrations of nickel; in fact the current values increased 
(see curves 6-8, Fig. 1). Extrapolation of the straight line portion of the 
curves to Y-axis indicated that initial values of the limiting current should 
be much greater than the observed values. 


The following Ilkovic’s equation gives the theoretical values of limit- 
ing current ig for different concentration of an electro-active substance: 


ig = 605 


where C is the concentration of the substance in millimolars; n, the 
number of electrons involved in the electrode reaction; D, the diffusion 
coefficient of the ion in cm.” sec.-!; m, the rate of flow of mercury in mgm./ 
sec.; and ¢, the time in secs. for 1 drop of mercury flowing through the 
capillary. Table II gives the observed values and the values of ig calculated 


from the above equation for different concentrations of Ni‘*; in these, 
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Fic. 2. Variation of the limiting current with the concentration of Ni+*. 
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D has been taken as 0-69 x 10-°cm.2/sec." The data in Table Il show 
clearly that the ig values do not follow Ilkovic’s equation above 15 mM. 
concentration (see also Fig. 2). It is interesting to note from these data that 


TABLE II 


Applicability of Ilkovic’s equation to the limiting current of Ni** 


ig calculated 


Concentration of from Ilkovic’s ig observed Deviation 
Ni*+* in mM. equation in in » amps. per cent. 
amps. 
(1) (2) (3) (4) 
2°5 7-67 7:7 1-0 
5-0 15-35 15-3 0:3 
7°5 20-02 24-1 4-6 
10-0 30-7 31-8 3°5 
15-0 46-05 44-3 3°8 
20-0 61-4 51-5 16-1 
(61-0) 
25-0 76°75 54-9 29-8 
(75-0) 
30-0 92-1 57-4 37-6 
(90-0) 
40-0 122-8 61-1 50-2 
(117-0) 


The values of i, given in the paranthesis refer to those obtained by extrapolation of 
straight line portion of amperometric curves to the Y-axis (see Fig. 1). 


the peculiar behaviour of the amperometric titration curves, viz., the non- 
linear variation of the current in initial stages of the titration, is closely 
associated with the non-applicability of Ilkovic’s equation at higher con- 
centrations. Since, however, the end points even at these high concentra- 
tions were very sharp, Ni** could be estimated polarographically upto 
40 mM. concentration using K,Fe (CN), as the titrant. 


j 
Aq 
= 
Pit 


96 H. C. SARASWAT AND S. K. D. AGARWAL 


SUMMARY 


The composition of nickel ferrocyanide was investigated following the 
amperometric titration of Ni*+ up to 40 mM. concentration against potassium 
ferrocyanide; it appeared to be 3 NiK, [Fe (CN),], Nis [Fe (CN),] or 
Ni;K, [Fe (CN),],. The data showed the possibility of estimating Ni+ 
polarographically within an experimental error of 1 per cent. 
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Received December 14, 1956 
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CARAMEL is a mixture of nitrogen-free organic compounds and is one of 
the colouring matters that are produced during the manufacture of cane 
sugar, chiefly due to overheating of sugars.'. Caramel not only affects the 
colour of sugar but also retards crystallisation.” It is manufactured on a 
large scale for purposes of colouring liquors, sweetmeats and gravies.*+4 


There are three main components in cane juice which are responsible 
for the production of caramel during the manufacture of cane sugar: these 
are sucrose, glucose and fructose. The process of conversion of sucrose 
into caramel acquires appreciable velocity only at pH S 9,° while in the 
manufacture of cane sugar, the different processes are operated, after the 
initial purification of the juice, nearly at neutral pH. Cane juice contains 
0-5 to 2 percent. of a mixture of glucose and fructose; these have been 
found to be the most important sources for the production of the colouring 
matter during sugar manufacture. The production of caramel from glu- 
cose has been studied by Geils,’ while the same from fructose has been exam- 
ined by Ripp® and others. 


The exact chemical composition of caramel is not yet available. 
Schumaker and Buchanan® reported that when solid sucrose was heated at 
about 200° C., three periods of foaming occurred, indicating the splitting 
of water at each stage and resulting in the formation of anhydrides with the 
empirical formule, and C3gH;9O.;. Recently, 
Luckesch!® chromatographically analysed caramel prepared from sucrose and 
reported that it contained glucose, fructose, glyceryl aldehyde, dihydroxy 
acetone, pyruvic acid, hydroxy methyl furfural and methyl glyoxal. Authors’ 
chromatographic studies" of the caramels prepared by heating glucose, 
fructose, maltose, sucrose and liquid glucose, in alkaline solutions, indicated 
the presence of four reducing compounds constituting caramel, with the 
following Ry (x 100) values: 0, 1-4, 4:4 and 6-5 (standard reference value 
= levulose — 13-6). The observation of a spot on the reference line on 
the chromatogram (Ry value = 0), which could not be separated from the 
base line even by running the chromatogram for more than 5 days, indicated 
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the presence of one or more reducing compounds of very large molecular 
weight." Authors’ chromatographic analysis indicated the presence of no 
acids in the caramels. 


Although attempts as mentioned above have been made to elucidate 
the composition of caramel, very few investigations seem to have been 
undertaken to study the kinetics or the mechanism of the production of 
caramel. The present communication reports data on chromatographic 
studies on the kinetics of the production of caramel from glucose. 


§2. EARLIER WORK 


The first attempt to investigate the kinetics of the development of 
colour in sugar solutions seems to have been made by Doss and Ghosh.!® 
These authors employed a spectro-photometer (Spekker) and studied the 
time development of colour in the invert sugar in alkaline solutions; they 
observed that “‘ double the concentration of invert sugar in the reaction 
mixture nearly doubled the rate of colour development”; and inferred, 
therefore, that the colour development in alkaline solutions at ordinary 
temperatures is essentially of the first order. The present authors'® also 
studied spectrophotometrically (using Unicam §.P.350 D.G. spectrophoto- 
meter) the development of colour in alkaline glucose solutions, at different 
temperatures in the range 35°-70°C.; and found that the following 
equation was applicable to the data on the time development of the colour 
in glucose at temperatures greater than 50°C.: 


a—r=ae*t 


where a is the total initial concentration of the substance undergoing change 
or the value of the optical density proportional to the concentration thereof; 
r, the amount decomposed upto time t; and k, the characteristic velocity 
constant of a first order reaction. This observation suggested that the pro- 
duction of colour or caramel from glucose, when heated in alkaiine solutions 
at temperatures above 50°C., is first order in nature. Employing the 
familiar Arhenius equation, the apparent energy of activation of the process 
involved has been computed from the data obtained at different temperatures ; 
an average value of this appeared to be 11-9 kcal. 


An interesting observation recorded by the present authors!*® and also 
by Doss and Ghosh,!? was that the initial portion of the curves represent- 
ing the variation of the optical density with time was concave in shape towards 
the X-axis (see Figs. 1 and 2); this feature is clearly seen from the data 
obtained at low temperatures, say 40° C. (see Curve 1, Fig. 2). At higher 
temperatures the concavity is noticeable only at very low time intervals. 
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§ 3. PRESENT INVESTIGATION 


The data in Figs. | and 2, especially the observed concavity towards 
X-axis of the optical density — time curves, suggested the following two 
important deductions’: 
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Fic. 1. Variation of the optical density with time during caramelisation of invert sugar 
(cf. Ref. 12). 


Fic. 2. Variation of the optical density with time during caramelisation of glucose at differ- 
ent temperatures (cf. Ref. 13). 


(i) the substance or one of the substances produced initially aids or 
catalyses the decomposition of glucose into different substances composing 
the caramel; or/and 


(ii) the substance or substances produced initially undergo further 
reaction presumably obeying a first order law, to give the final products; 
that is, the conversion of glucose into caramel is apparently through the 
formation of intermediate compound(s). 

These deductions have now been examined, in some detail, employing 
the filter paper chromatographic technique. The procedure involved the 
chromatographic investigation of the reaction mixture at different intervals 
of time during the progress of the reaction. 


§4. EXPERIMENTAL 


Analar samples of glucose and sodium hydroxide were used. Known 
quantities of solutions of these substances, of desired strength, were mixed 
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together such that the final solution was of a pH value in the range 10-13, 
About 100c.c. of the mixed solution was taken into a glass vessel which 
was immersed in a paraffiin oil thermostat. The temperature of this last 
could be raised to and maintained constant (+: 0-1°) at any desired value. 


During the development of colour in the alkaline glucose solutions at 
a fixed temperature, a small quantity of the reaction mixture was withdrawn 
at regular intervals of time and was used for chromatographic analysis. It 
is to be pointed out that the portion of the reaction mixture withdrawn for 
analysis was quickly cooled to the room temperature and thus the progress 
of the reaction in the sample got practically quenched. 


The method employed for chromatographic analysis was essentially the 
same as used by Patridge." Strips, of suitable width, of Whatman filter 
paper No.1, were used. The reaction mixture under examination was 
spotted on a (reference) line drawn parallel to one edge of the filter paper 
strip. The solvent used was the top layer of a mixture of butanol, ethanol 
and water (50: 10:40 by volume). The bottom layer was used for saturat- 
ing the chamber in which the filter paper strips were suspended. After the 
chromatograms were run for 70-100 hrs., they were dried and washed with 
a solution of silver nitrate in acetone. They were later dried again and 
sprayed with alcoholic KOH solution, when brown or dark brown spots 
referring to the reducing compounds present in the reaction mixture, became 
apparent. The chromatograms were later treated with ammonia, washed 
with continuous flow of water for 2 hrs. and finally dried. 


§5. RESULTS 


Figures 3 and 4 refer to two typical chromatograms. These refer to 
reaction mixtures drawn at regular intervals from the reaction vessel in which 


TIME Wm MINUTES TIME IN MINUTES 
° 30 o t 2 5 10 20 30 


Fic. 3. Fic. 4. 
Fic. 3. Chromatographic studies of the composition of the reaction mixture at different time 
intervals during caramelisation of glucose at 40°C. 
Fic. 4. Chromatographic studies of the composition of the reaction mixture at different 
time intervals during caramelisation of glucose at 70° C, 
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glucose has been heated with alkali solution at 40° and 70° C. respectively. 
The running time for the chromatograms of Fig. 4 was 72 hrs., while the 
same for the others given in Fig. 3 was 90 hrs. The Ry values for the spots 
observed in Figs. 3 and 4 are returned in Table 1; these were computed by 
comparison of the spots developed for known sugars in a previous run on 


a separate strip of filter paper. 


TABLE | 
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Ry values of the spots observed in the chromatographic studies on 
the kinetics of caramel formation 


Time in Minutes 


Spot No. 
0 30 100 200 
Temperature of Reaction Vessel = 40° C. 
| 0 0 0 
2 2°1 2°8 
3 4-5 4-5 4-3 
4 6°5 
5 9°5 9-5 9-5 
6 12-8 12-9 13-0 
7 14-9 15-0 15-0 
0 1 2 5 10 20 30 
Temperature of Reaction Vessel = 70° C. 
1 0 0 0 0 
2 2°8 2°8 
3 4-3 4-4 4-4 4-4 
4 7°8 7°8 7°8 7°8 7°8 7°8 
5 9°5 9-5 9-5 9-5 9-5 
6 12-8 12-8 12-2 12-2 12:2 
7 14-2 14-2 14-2 14-0 
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The chromatograms taken with the reaction mixture obtained in the 
beginning of the reaction, i.e., at zero interval of time, indicated the presence 
of only one spot corresponding to glucose [with an Ry (x 100) value of 9-5], 
The chromatograms referring to the reaction mixture drawn at increasing 
intervals of time, showed spots other than that of glucose. It was interesting 
to note that the reaction mixture obtained at 2 mts. intervals of time contain- 
ed, as indicated by chromatographic analysis (Fig. 4), three compounds, 
excluding glucose, the Ry (x 100) values of these are approximately as 
follows: 3-0, 4-1 and 12. The production of a compound with an Ry; 
value of 12 referring to fructose (Ry value 13 was particularly instructive 
in light of the finding due to Lobry de Bruyn and van Ekenstein® who showed 
that during heating of sugars in alkaline solutions for caramelisation, 
equilibrium mixture of glucose and fructose (as also mannose) are produced; 
in accord with this is the observation of the authors!! that whatever may be the 
sugar employed initially for caramelisation, the same reducing compounds 
have been detected in chromatographic analysis of the caramels. Further, 
the data in Figs. 3 and 4 show that the spaces where the original mixture was 
placed on the reference line for chromatographic analysis, acquired colour 
only after 25 mts. in the reaction mixture maintained at 40° C. (Fig. 3) and 
2 mts. interval of time in the reaction mixture maintained at 70° C. (Fig. 4). 
This spot (Ry value = 0) on the reference line, which could not be separated 
from this last even by running the chromatogram for 5 or 6 days indicated, 
as mentioned above (§1), the presence of reducing substance(s) of high 
molecular weight. This observation suggests that the substances compos- 


ing caramel are produced after the initial generation of certain substances 
with low molecular weight. 


§6. DiIsCuSSION 


The following is the resume of the important findings on the kinetics 
of the production of caramel: 


(i) During heating of glucose and other sugars in alkaline solutions, 
compounds of low molecular weight are produced initially, 


(ji) Lobry de Bruyn and van Ekenstein transformation takes place; 
(iii) the production of caramel involves autocatalysis; and 


(iv) the reaction involves a step which controls the time development 
of colour and which obeys a first order law. 


These observations lead the authors to suggest the following mechanisms: 


(A) Heating of glucose or other sugars in alkaline olutions, leads to 
the initial production of an equilibrium mixture of glucose, fructose and 
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mannose; these apparently condense together to form an _ intermediate 
compound of very high molecular weight, which decomposes unimolecu- 
larly to yield the compounds [with Ry (x 100) values 0, 1-3, 4-4 and 6:6, 
in reference to levulose — 13-6] constituting caramel, one of which acts as 
a catalyst to the establishment of Lobry de Bruyn and van Ekenstein 
transformation. 


(B) Alternatively, the rate determining step in the process is likely to 
be the establishment of Lobry de Bruyn and van Ekenstein transformation 
and the autocatalytic nature is related to the decomposition of the conden- 
sation product contemplated in (A). 


The observation of the autocatalytic nature of the reaction upto 30 mts. 
at 40°C. and upto 2 mts. at 70°C. suggests that the mechanism (A) is to 
be preferred. Conclusive evidence is, however, lacking. 


§ 7. SUMMARY 


The kinetics of production of caramel which is a mixture of nitrogen- 
free organic compounds and which offers colour to sugar crystals and 
retards the rate of crystallisation, have been investigated by employing filter 
paper chromatographic technique. Chromatograms of the reaction mixture, 
viz., glucose heated in alkaline solutions, at different intervals of time during 
the progress of the reaction showed the initial production of compounds 
of low molecular weight, which condensed to form an intermediate compound 
of very high molecular weight; this last appeared to decompose unimole- 
cularly to give the products composing caramel. 
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NEW SOLUTIONS FOR FINITE DEFORMATION 


By B. R. Sern, F.A.Sc. 
(Kharagpur) 


Received January 10, 1957 
|. INTRODUCTION 


IN recent years the theory of finite deformation has received fresh impetus 
from a number of exact solutions that have been obtained for imcompressi- 
ble elastic bodies without assuming any form of the strain energy function. 
The exact solutions obtained for compressible bodies have not been very 
large, and in general they have been of the form obtained by B. R. Seth. 
Extensive references to existing literature on the subject are given in the 


Presidential Address to the Section of Mathematics at the 42nd Session 
of the Indian Science Congress.! 


Very few attempts have been made to obtain solutions for the finite 
bending of plates into shells. The earliest exact solution on the basis of a 
linear stress strain tensor law was obtained by Seth? for a rectangular plate 
bent into a right cylindrical shell. The form of this solution has been used 
by many workers. In the present paper an attempt has been made to dis- 
cuss solutions which may be used to discuss the finite bending of plates 
into spherical shells. A linear stress strain tensor law is again assumed, 
and this gives sufficiently good results for technical applications. 


In cylindrical co-ordinates it is found that, if we take 
u=r—¢(r,z), v=9—AI, w=z—f(r? 4+ 2%), 


all the equations of finite deformation and boundary conditions on a 
spherical shell can be satisfied only if 


d=ar, f(r?+2)=— (r? + 


where a, A, c, are known constants and c = (1 — 2o)/(1 — a) is an elastic 
constant. 
2. Frnire STRAIN COMPONENTS IN CYLINDRICAL CO-ORDINATES 


For all applications finite strain components should be referred to the 
strained state of the body. Let *x be the unstrained and x the strained 
co-ordinates of the deformed body with reference to a fixed frame with 
jg and gj; as metric tensors. The strain tensor ej; is given by? 

= Sig — ap (2.1) 
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Putting ut = xt —+x, (2.1) becomes 
= Bij — + + — uP, (2.2) 
In cylindrical co-ordinates 
ix = 19, 0, =r, 90, z, 


r¢ @ @ @ 
=] ro? 9 r? Of. (2.3) 
0 0 1 0 0 l 


The tensor components «; can be replaced by the physical components 
e4; by the relation* 


= (gitgii)te;;. (2.4) 


Thus we get the components of finite strain in cylindrical co-ordinates as® 
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1 dw ny) dv 


1 du du dw dv 1 dw dw 


~ 450 dz r dz’ 


ow du du 2 OV ov ow dW 
ES + | 
dv OV OV dv 
1 du dv 


du du dv ov 1 dw dw 


We have now to use the semi-inverse method of assuming a possible 
set of displacements. For the type of problem in view we should take 
u=r— o(r,z),v = 9— Ad,w =z — f(r? + 2%), (2.6) 


where f and ¢ are arbitrary functions to be determined and A is a constant. 
The corresponding strain components are found to be 


Crp = £(1 — — 2r2f’?, = Ky — 34 2° 


Cog = 4(1 — — 42°"), = — (Grbz + 4r2f"), 
Coz = 0, Crg = 0, (2.7) 
where 


K,=34—(1—A), K,=1—2(1 — A) — (2.8) 


by =F, R? = r? + 2?. 
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3. STRESS-STRAIN RELATIONS AND BODy-STRESS EQUATIONS 
We assume that 
Tig = + (3.1) 


The body-stress equations of equilibrium are 


Ti, j = 9, 
and 


4=i—K,—2¢* + 2) f"?—4K,% — (Gr? + 2"). (3.3) 


Substituting the value of e;; from (2.7) we see that the first equation gives 


0 19 9 9 
yp LAS — J — J — + 
2f > (¢ 1 
(1 2Ky + Ky — pd,” + | = 0. (3.4) 
The second is identically satisfied and the third gives 
Ad — J — 8u J — (by? + 


2 bz bz 
+ [2 r =0 (3.5) 
These will be satisfied by 


P = Ad — f — 84 f — (by? + 


=a constant, R? = r? + 2?, (3.6) 
provided 
d (bz 
or 
‘ 
aC — 2Ky + )- br? E (3.8) 
From (3.7) we get 
— (3.9) 


which gives 
= F [x (z) + logr], (3.10) 


wh 
Sul 


Fr 
Ce 
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where x and F are arbitrary functions of z and x (z) + logr respectively. 


Substituting this value of @ in (3.8) we get 


(1 — 2Ky) r? + K,F? — r?F’?y’” (z) — F’? = 0. 


From (3.6) we see that 
2 


WK, + (ZA + (Gp? + = M (r? + 2°), 


(3.11) 


(3.12) 


M being an arbitrary function of (r* + =*). Substituting from (3.10) we get 


Comparing (3.11) with (3.13) we get 


M = — 3A(1 — 2K,) 


A 


x" = 
Thus 

x = — log (cz + D), 
and 

[log p| 
Also 


r2c2 
F’? = (1 — 2K,) r? + K,F%, 

which shows that c= 0 and 

F = aeér, omitting a constant, 

= af, 
where 
1 — 2Ky + K,a? — a? = 0, 

or 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


| : 
a 
) 
a 
| 
=1l, 
A=1, 
) Thus we see that the value of ¢ is 
- 
1—A 
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4. DETERMINATION OF f 
Substituting the value of 4 from (3.18) in (3.6) we get, 
A[l — Ko — — —2R*f"?] — — 8p 
— 8u R2f’f"dR* = a constant. (4.1) 
By differentiating (4.1) with respect to X = R? we get 


20) A+ = 0, (4.2) 
which gives 
f= mre, (4.4) 


neglecting a constant. 
The components of displacements now take the form 


+ (4.5) 


u=r—ar, v = 0— A, 


where c is the elastic constant (1 — 20)/(1 — a) and a, A, c, are constants 
to be determined from the boundary conditions and equation (3.17). 
The non-vanishing stress components are 


Trp = AD + 
A(L — Ko — $Kya? — + (1 — a2) 
2Ac,2R-? 
Tyg = A(L — Ko — $Kya? — 4a”) + w (Ky — $Ka?) — 2Ac,*R-, 
Tzz = A(1 — Ko — — + — 2Ac?R-* 
4uz*c,*R-? | 
Tre = — 4urzey2R-? 
Tre = 0, T oz = 0. (4.6) 
5. BOUNDARY CONDITIONS FOR A SPHERICAL SHELL 


If N denotes the direction of the normal at any point of the shell, the 
boundary tractions in the direction of r, 9, z corresponding to the general 
set of displacements given in (2.6) are 


ou 


a 
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— (rb + (5.1) 

— + 22), | (5.2) 
N, = 0. (5.3) 


substituting the values of ¢ and f’* from (3.18) and (4.3), we get 
RNy = r [A(] — Ky — — + — pa? 


— 2(A + 2p) (5.4) 
RNz = z[A(1 — Ko — — + 
2 (A + 2p) (5.5) 


If we assume that the inner surface R = a of the shell is acted upon by 
a normal traction of the amount }ya*, we get 


— Ko — — $a*) + p — 2(A + 2p) = (5.6) 


Assuming the external surface R = b to be acted on by a traction of amount 
po?z/b in the direction of the axis, we get 


A(1 — Ky — $K,02 — + —2(A4+ c,2b-% = pa? (5.7) 


These give 


and 
Ja? | t 2-0] = (1-0). (5.9) 
Also from (3.17) we have 
K,=1, (5.10) 
1—A 


(ii) a2 = 


A’ 


2) 
a 
3) 
7 
| > 


112 B. R. SETH 


2 (px — ar 
(5.11) 
(ab)*° 


In (ii) A may be determined from (5.9) and substituted in (5.8) to get c,?. 


Assuming that the plate is bent into a spherical shell with a plane rim, the 
tractions on the rim are 


Zy = Tez = — (5.12) 
Zz = Tzz = A(l — Ko — 4Kya? — 4a?) + — 
— 4uc,2z*R-? ©), (5.13) 
Their statical resultant per unit length of the rim z = Zp, is 


rdr. = 4A(1 — Ky — — 402) (6? — a2) 


Ac,” 


“gece — gi) + a’). 


This is the force which must be applied on the rim to keep the plate in the 
spherical form. 


6. SUMMARY 


The use of finite components of strain shows that a plate can be bent 
into a spherical shell if suitable tractions are applied to the inner and outer 


surfaces and to the plane rim. The corresponding components of displace- 
ment have been obtained. 
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Received December 26, 1956 
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INTRODUCTION 


EsTIMATION of copper in solutions!,? and in a number of alloys** and 
metals,®.* as also in semi-conductors’ has been the subject of recent interest. 
For this purpose a series of reagents like dithizone,® thiobenzamide,' 2-methyl- 
8-quinolinol,® nitrilotriacetic acid,’ diethyldithiocarbamate,® ethylenedi- 
aminetetra-acetic acid (EDTA),":!* etc., have been developed for colorimetric 
estimation of copper. A resumé on the work done, during a period of two 
years (from August 1953 to August 1955) amounting to a few scores of 
papers, has been given by Mellon and Boltz.'* 


It is known that EDTA forms a soluble complex with copper"; using 
this property, simple methods for estimation of copper in steel, bronze and 
other alloys, have been developed. Thus Sedivec and Vasak” determined 
copper in steel by complexing the mixture with EDTA and later precipitat- 
ing copper with diethyldithiocarbamate. Similarly Byrn and Robertson!® 
precipitated copper quantitatively in the presence of EDTA from a solution 
of nickel copper alloy or bronze by 2-(0-hydroxy phenyl) Benzole. 
Direct spectrophotometric methods for determining copper with EDTA 
have also been developed. Malmstadt and Gohrbandt'’? recorded the 
absorption of copper nitrate versenate and copper versenate solutions in 
the ultraviolet region and used the absorption at A = 320 mp for estima- 
tion of copper with versenate and vice-versa. Similarly in the near infra red 
region Hughes and Mortell'® showed that copper versenate has its maximum 
absorption between A = 700-750 mu; employing the absorption of copper 
at A= 745 mp, Underwood” titrated copper solution with EDTA; and 
Nielsch and Boltz!* estimated copper with EDTA using the absorption 
maximum at A = 720 mp. 


In the course of a series of investigations on the easy estimation of 
calcium in sugar solutions and sugarcane juices, which is a problem in sugar 
technological laboratories, a new method for the estimation of copper 
spectrophotometrically has been evolved, This method appeared usefuy 
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on account of the possibility of using the characteristic absorption of 
copper in the visible region (at A = 630 mp, vide infra). It may be added 
that there have been similar simple methods worked out (vide supra) which 
involved the use of absorption of copper in the ultraviolet region or in the 
near infrared. Non-availability of the ultraviolet and other spectrophoto- 
meters in our Laboratory set a limitation to the use of these methods. The 
present communication describes a method for estimation of copper in 
ammoniacal medium by disodium salt of ethylenediaminetetra-acetic acid 
using Unicam SP 350 DG Spectrophotometer, which has a _ continuous 
variation of wavelength in the region 400-700 mp. 


EXPERIMENTAL 


The spectrophotometer employed, viz., Unicam SP 350 DG Spectro- 
photometer has an one cm. light path. 


The substances used were of the following grade: Copper sulphate 
used was of BDH Analar quality; this was purified by recrystallisation in 
distilled water. 


The sodium versenate employed was of Merck quality. 


The method for carrying out the spectrophotometric titrations was 
as follows: Solutions were made in such a way that the final volume (25 ml.) 
contained a fixed quantity of copper and ammonia and varying quantities 
of the titrant, viz., the sodium versenate. The absorption of these solutions 
was studied in a few cases, in the range 400-700 my and in the subsequent 
experiments at A = 630 mp (vide infra). 

RESULTS AND DISCUSSION 

Figure | gives a typical series of results on the absorption of ammoniacal 
copper solution in the range 400-700 mp. In this, curves 1, 2 and 3 refer 
respectively to 20, 16 and 10 mM. concentration of copper solution con- 
taining 1-15 M. of ammonia. It was interesting to note that ammoniacal 
copper solution exhibits well-defined characteristic maximum at A = 630 mp; 
the absorption at this wavelength obeys Beer’s law (see Fig. 2) with a value 
of 50 mole litre cm. for the corresponding constant. Further, it was 
instructive to note that the absorption of copper versenate at A = 630 mp 
also obeyed the Beer’s law (Fig. 3); the Beer’s constant for copper versenate 
under the experimental condition, viz., in the ammoniacal medium was 16:7 
moles! litre cm.! This appreciable difference between the  Beer’s 
constants of copper and copper versenate solutions facilitated the execution 
of the titrations under investigation. 

Figure. 4 gives a typical set of curves representing the variation of the 
absorbency of ammoniacal copper at A = 630 mp during the progress of the 
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Fic. 4. Spectrometric titration of Cu+t+ with sodium versenate at 4 = 630 mu. 
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titration. It was instructive to note that optical density of copper decreased 
progressively with the addition of versenate till the end point; further 
addition of versenate effected little the optical density. It may be noted that 
the end points were very sharp; the error involved in the determination of 
copper by versenate is less than 1%. 


SUMMARY 
Employing the absorption at A = 630 mu of copper ammonium com- 


plex, a new method for determination of Cu*+ by sodium versenate was 

described. The experimental error involved in the method was less than 1%. 
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1. INTRODUCTION 


RECENTLY Van Niekerk and Schoening! have studied the structures of a 
number of metal acetates. They have stated that direct metal bonds exist 
in some of these compounds. Magnesium acetate tetrahydrate belongs 
to the same group as those studied by Niekerk et al., and the present paper 
gives a report of the investigation of its structure by X-rays. 


2. EXPERIMENTAL 


The crystals of magnesium acetate were obtained from its aqueous 
solution by slow evaporation. The crystals were needle-shaped, the long 
axis being the c-axis. The unit cell dimensions, as determined from rotation 
and Weissenberg photographs, are 


a=4-75A, b=11-799A, c=8-52A and B= 94° 54’. 


The axial ratios a:b: c = 0-402 : 1: 0-719 are in good agreement with the 
previously reported morphological values a:b: c = 0-403: 1 :0-7128 and 
B = 95° 37’ (Groth?). Taking the density to be 1-453 (Hand-book of 
Chemistry and Physics) the number of molecules per unit cell is 1-98, i.e., 2. 
There were no systematic absences in the hk/ reflections showing that the 
unit cell is primitive. The list of reflecting planes reveals that 0 k 0 is absent 
when k is odd, and h0O/ when / is odd. The space-group is therefore 
Cy,° — P2,/c. 


Using the multiple film technique, the zero layer Weissenberg photo- 
graphs were taken on suitable crystal specimens rotating about the a- and b- 
axes using Cu-K, radiation. The intensities of the reflections were estimated 
visually and corrected in the usual way for the Lorentz and polarisation 
factors. In view of the small size of the crystals used no corrections were 
considered necessary for absorption and extinction. 
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3. DETERMINATION OF THE STRUCTURE 


Since the general point position in the P2,/c is four-fold and as the unit 
cell contains only two molecules, it follows that the magnesium atoms must 
be on symmetry centres. The magnesium position was therefore taken at 
(0, 0, 0). Since the cell dimensions of magnesium acetate and nickel acetate 
are almost identical, and since they belong to the space-group C,;° with the 
metal atoms at special positions, the signs of the structure factors of mag- 
nesium acetate were deduced from those for nickel acctate, taking into 
consideration the contribution due to the metal atoms. With the signs 
so obtained a trial Fourier projection on the bc plane was made. The 
projection thus obtained clearly showed the positions of all the atoms and 
the peaks were well resolved. The maxima were located for all the atoms 
by Booth’s method,* and the y and z co-ordinates of the atoms were estimat- 
ed. The process was repeated for the ac plane, and a Fourier projection was 
made. Figs. | and 2 give the Fourier projections on bc and ac planes res- 
pectively. From these projections the co-ordinates were estimated (Table J), 


Fic. 1. Fourier projection of the electron density on (ac) plane. 
With these co-ordinates, structure factors were calculated by allowing for 
the temperature corrections, The f-values for carbon were taken from 
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H,0(2) + 


Fic. 2. Fourier projection of the electron density on (bc) plane. 


TABLE 


Atomic parameters 


0-000 0-000 0-000 
—0-250 —0-150 —0-015 
0-067 —0-267 0-150 
—0-283 0-083 0-133 
0-200 —0-050 0-217 
—0-183 —0-250 0-083 


0-300 0-367 0-001 


Robertson,! for oxygen they were proportionately deduced from Robert- 
sons’s values, and for magnesium atom an isotropic temperature factor with 
B= 1-2 was used. In Table II, the calculated and observed structure 
factor values are given. The agreement is good and the Reliability Factor 
R comes out to be 0-13 for [((0k J] zone and 0-21 for [(h0/)] zone. The 
bond length and the angles are given in Table III, which were deduced from 
the atomic co-ordinates given in Table I. 
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TABLE II 
| | | | 

hei] Fo | & F, Fo | & & 2 F, Fy 

10 +45 3 0451{ +4 4/023! 4 
20 0] +3 30;046)/! +2 6 | 0 12 4 1 3 
30 417 | 17! ll | 011265 +83 4 
40 +3@ 2% 0 5 2! +09 | 0 1814 +7 5 
5 0 0 +9 5 0 53 | +1 19 0 13 5 + 4 4 
o2 0/ -7 8/1 02 +30 | 2% 
04 0} +17 17) 0 413 13 | 1 0 4 | +22 | 19 
«2 5/0 656] +4 6/106 +5 8 
08 0} +21 46 42 4/1 08 +19 | 12 
00 4/ +2 6/061 +13 | 211 5 
00 6) +32 2% 62 +8 2/1 06 
00 +15 +65 6,108); +2 4 
0010 +17 12 0 64 43) 02 -7 7 
1; +387 3 66! +12 9 204) +4/ 2B 
01 +24 2 68!+4+9/] 8 2 O06 | +12 | 15 
01 4] +3 3/2 0232) +20 | 16 
01 6/ +5 4'075 +16 | 11 2 08 +9 7 
ot t | #8 10 0 76 +4 4 0 2 7 
01 8/+5 40] +71] 9 
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TaBLe []—(Contd.) 


Fo | & 2 Fe | Fo hkl F, Fo 
| 

01 +9 6,078] -2 4 3 06) +28 | 18 
110 -2 4 0 6 3 08) +8 4 
111 +3 4 O 81 +14 | 14 3 0 2 | +34 | 2% 
2 3 -10 8 0 83! +38 2 83 06! +8 8 
24: 20 0 8 +11 9 +16 | 
2 +11 | 12 85] -5 9 402)| +64 5 
2 6, +10 12 8 4041 +9 7 
27\|+4 2/0 88s +8 5 4 06. +15 8 
2 8) +10 6 +4 8 408 +2 3 
3041 +17 15 | 0 9 | 16 | 4 0 2) +410 7 
3 2 -% | 983! +18 4 04! +12 6 
+e! «lie oe | «le 
3 4 2/0 96 +6 8 4 08 | +8 8 
3 5 +4 + +6 7 5 5 
3 6 +16 +121 4 5 041] 411 8 
3 8 | -12 10/0102 +6 -83 4 
3 +8 7/0105 +1 2 5 -6 5 
310, +7 8 | +413 9 5 O02] +8 6 
4 1 +9 -8, 10 5 04|+6) 6 
42, +8 5|/oh 4 -5 8 5 06 7 5 
43, -15 | 5 6 
44: -3 5/0122 +9 UW | 
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TABLE III 
Bond lengths and angles 


O,-O, .. 2:23A Mg-O, 2-1A 


O0,-C, ee oe 33 Mg-H,0O,,) 
0,-C, 


C,-C, 


2-08 
Mg-H,O,,, 


O,-Mg-H,0,, 90° 48’ 


88° 48’ 


4. CONCLUSIONS 


The structure is anologous to that of nickel acetate. The two acetate 
groups are planar and are distributed with the magnesium atom at a 
centre of symmetry. Each magnesium atom is surrounded by four water 
molecules H,O,1,, HO;2), H2O,2), and oxygen atoms O, and O, (bars 
denoting the atoms belonging to other acetate group, not shown in the 
figure). The remaining oxygen atom Oj.) is linked to H,O,.). The acetate 
group belonging to the magnesium atom situated at (0, 4, 4) is also linked 


to the acetate group at (0, 0, 0) by hydrogen bond, as in the case of nickel 
acetate. 


The bond distances between the magnesium atom and water molecules 
are 2-08 and 2-07 and the distance Mg-O, is 2-11 A. These are in good 
agreement with the values obtained in the case of nickel acetate. The 
O, - C,, and O, - C, are found to be 1-33 and 1-25, while in the case of nickel 
acetate they are reported to be nearly equal. 


5. SUMMARY 


The crystal structure of magnesium acetate has been completely 
determined. The crystals are monoclinic with space-group C,,° — P2,/c 
and two molecules in a unit cell with edges a = 4:75 A, b = 11-79 A, c= 
8-52 A and 8 = 94° 54’. The positions of the atoms were determined from 
electron density projections on the bc and ac planes. From these the various 
bond lengths and angles were calculated. 
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l. 


INTRODUCTION 


UNLIKE other components of the connective tissue, relatively few structural 
studies have been reported on elastin. The earliest to investigate the X-ray 
diffraction pattern of elastin was Kolpak (1935) who found that the unstretch- 
ed ligament yielded diffuse rings, but on stretching these gave a typical collagen 
pattern with meridional ares corresponding to a spacing of 3 A and equatorial 
spots corresponding to spacings of 11-5 A, 5-9 A and 4-6 A. Astbury (1938, 
1940) attributed these results to the presence of small amounts of collagen 
as impurity in the specimens. According to him, the unpurified fibres of 
ligamentum nuche give only an amorphous pattern, but this may be super- 
imposed with the pattern of imperfectly oriented collagen fibres. On stretch- 
ing, the collagen pattern becomes well oriented. On the other hand, he 
found that the collagen pattern disappeared if the fibres were treated with 
hot water over a steam-bath for about a week. On stretching, the amor- 
phous pattern of elastin also showed some signs of orientation parallel to the 
direction of stretching. Astbury suggested that elastin might be a member 
of the collagen group whose thermal transformation temperature is below 
room temperature, but no conclusive evidence was available at that time to 
decide the question. Bear (1944) recorded small angle patterns of collagen 
faintly with beef ligament and ascribed these to the presence of very small 
amounts of collagen as impurity. He suggested that elastin must be ex- 
cluded from the collagen family (Bear, 1952). Electron microscopic studies 
(Wolpers, 1944; Gross 1949, 1951; Franchi and de Robertis, 1951) also 
did not reveal the 640 A spacing typical of collagen or any other fine structure 
in elastin. The present situation regarding the structure of elastin, as 
revealed by X-ray and other physical methods, may be summarised as 
follows (Kendrew, 1953): ‘‘ Indeed the two main groups of the Astbury 
classification appear to comprehend virtually all known structure proteins 
whose X-ray pattern has been examined, with the exception of elastin (a 
component of some types of elastic tissue) and actin in muscle.” 


Among the other chemical and other studies on elastin, the following 
may be mentioned in particular. Lloyd and Garrod (1948) compared the 
thermoelastic behaviour of elastin and collagen and considered that the 
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high proportion of non-polar side chains present in elastin would account 
for the rubber-like properties. Hall eta/. (1952) suggested that the 
elastin system at least in parts is made up of linear aggregates of a corpus- 
cular protein unit of low molecular weight called proelastin, analogous to 
procollagen. Banga and Schuber (1953) have suggested that elastin is a 
glycoprotein, the fibrous character of which is determined by a carbohydrate 
group built into the molecule. Later Banga (1953) showed that the reducing 
substances split off during elastolysis could not be ascribed to carbohydrate, 
Partridge er al. (1955) have suggested that the substance of elementary elas- 
tin fibres is a chemically homogeneous protein, although it may exhibit 
heterogeneity at macromolecular level. The physical properties of the 
elastic tissue may be a result of the macromolecular structure of elastin 
comprising the elementary fibres and the way in which the fibres are cement- 
ed by collagenous structure and mucoproteins. 


During a recent series of studies in this laboratory on the effect of 
chemical reagents on collagen, it was observed that some of the treated speci- 
mens exhibited a very close similarity to elastin in their X-ray pattern as well 
as in the physical behaviour (unpublished). The change thus produced in 
collagen was found to be reversible, so that the treated specimen should be 
considered to be only a modified form of collagen. A careful investigation 
was therefore made of the relationship between collagen and elastin and the 
present paper is a report of these studies. It appears that elastin should 
definitely be classed with collagen and that its molecular structure is also 
basically built up on the triple chain structure of collagen (Ramachandran 
and Kartha, 1954, 1955, a, b; Ramachandran, 1956). 


2. X-RAY PATTERNS OF ELASTIN AND COLLAGEN 


The X-ray diffraction pattern of elastin supports the suggestion that 
elastin is closely related to collagen in its structure. The diffraction pattern 
of elastin, recorded in a 3cm. cylindrical camera is given in Fig.1. The 
pattern in Fig. 1 may be compared with the diffraction patterns of natural 
collagen fibres (Fig. 2), thermally shrunk collagen (Fig. 3) and gelatin 
(Fig. 4). The fibres of elastin were obtained from ligamentum nuche of cattle. 
The fibres were treated with 20 percent. sodium chloride solution, washed 
with distilled water and dehydrated with acetone. Further purification was 
effected by autoclaving with water. Collagen fibres were obtained from 
cattle achilles tendon from which soluble proteins were removed by extrac- 
tion with sodium chloride solution. B.D.H. gelatin was used for obtaining 
Fig.4. The thermally shrunk collagen specimen was prepared by heating 
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a collagen fibre well above the shrinkage temperature and then fixing it jn 
boiling alcohol. 


The prominent features of the diffraction pattern of elastin are two 
diffuse rings of spacings 4:4 A and 2:2 A. While the former was recorded 
by Astbury using a flat film, the latter has now been observed for the first 
time. The pattern recorded on a flat film is shown in Fig. 5 and it will be 
seen from this that there is in addition a central halo, with a faint maximum 


at about 9-10 A. This feature has also been commented upon by Astbury 
(1940). 


A comparison of Figs. | to 4 shows that the 2-2 A ring is common to 
all, although the other details differ. In fact, there is a close correspondence 
between the patterns of elastin and thermally shrunk collagen. However, 
in the patterns of collagen fibre and gelatin, several other reflections occur. 
Thus, the pattern of collagen fibre exhibits a sharp equatorial reflection at 
11 A to 12 A and meridional arcs at 2-86 A and 4 A, and two off-meridional 
spots corresponding to a spacing of about 7A. All these reflections occur 
in the form of rings in the pattern of gelatin. These features are not fully 
visible in Figs. 2 and 4since the central region has been overexposed to 
bring out the 2-2 A ring. Thus, gelatin at room temperature may be consi- 
dered to be a disoriented (polycrystalline) form of collagen. 


The pattern of thermally shrunk collagen differs significantly from 
those of both normal collagen and gelatin. While the occurrence of complete 
rings indicates the absence of orientation, as in gelatin, there is a clear 
indication that it is a more degraded form of collagen than gelatin, since the 


sharp rings in the pattern of the latter are all absent in the pattern of 
thermally shrunk collagen. 


The sharp reflections corresponding to definite spacings occur in the 
pattern of collagen fibres because of the high degree of order in its structure. 
This orde: is of two types: 


(a) the order arising from the arrangement of the triple-chain proto- 
fibrils in a lattice—the Jattice may be a simple hexagonal lattice, or a 
cylindrical lattice (Ramachandran and Sasisekharan, 1956) but the exact 
nature is not important for our discussion here. 


(b) the order arising from the helical arrangement of residues in an 
individual triple-chain protofibril. 


The former leads to the 11A to 12 A equatorial spacing, while the latter 
is responsible for the 2°86 A meridional arc and also the other spacings. 
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The 2:86 A corresponds to the tenth order of the repeat distance of 28-6 A 
spacing along the fibre axis. 


It is obvious that the two types of order must be interrelated. An 
examination of the way in which the triple-chain structure is built up from the 
amino-acid residues and is stabilised reveals that the exact value of 28-6A 
is not demanded by stereochemical or any other considerations connected 
with a single triple chain. A fairly wide variation, certainly from about 
98 A to31 A, is permissible. Thus, the occurrence of a sharp spacing of 
2:86 A must be attributed to the cross-links between different protofibrils, 
which would lead to the repeat spacing being the same in all the protofibrils, 
and thus producing a high degree of regularity. If for any reason these cross- 
links are ruptured, then one should observe two effects—(a) the equatorial 
reflection would no longer be sharp, because there is no reason why the 
protofibrils should be maintained at a constant distance apart, (b) the indivi- 
dual protofibrils may also have different values for the repeat spacing along 
the fibre axis, so that the 2-86 A arc would be rendered diffuse and may even 
disappear. This is precisely what happens as a result of thermal contraction. 
The sharp equatorial spots in collagen fibre, which correspond to the 11 A to 
12 Aring in cold gelatin, disappear on thermal contraction, and simultaneously 
the 2°86 A arc and the other sharp reflections or rings also disappear. 


The fact that the broad rings at 4-4 A and 2-2 A, which are observed 
with gelatin, are present in the diffraction pattern of the contracted collagen 
suggests that the essential triple helical structure is maintained, except that 
the lattice structure and regularity are lost in the latter. This suggestion is 
further supported by the fact that thermally contracted collagen may be 
readily brought to the normal state by stretching it at room temperature. 


As already mentioned, the diffraction pattern of elastin is closely similar 
to that of thermally contracted collagen, and so it is reasonable to suppose 
that the structure of elastin is also based on the triple-chain collagen structure, 
in which however the cross-linkages are not present. This happens because 
elastin contains very few residues with polar side-groups and in particular 
contains negligible amount of hydroxyproline (Partridge et al., 1955), which 
is the main cross-linking agent in collagen. 


3. THEORETICAL CONSIDERATIONS 


The above observations strongly suggest that the structure of elastin 
and collagen are closely analogous at the molecular level, and that they are 
both built up on the triple-helical arrangement proposed from this labora- 
tory (Ramachandran and Kartha, 1955; Ramachandran, 1956). This is 
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theoretically also very reasonable, for the amino-acid composition of elastin 
is similar to collagen in two important respects, although there are very 
few residues with polar side-chains in elastin. It contains slightly more than 
one-third of glycine and further the proline content is nearly the same as in 
collagen (Tristram, 1953). The triple-chain structure has the essential 
property that one-third of the a-carbon atoms cannot have any side-groups 
attached to them so that the glycine residues must form at least a third of 
the total number. This is characteristic both of collagen and elastin. So 
also the five-membered ring of the proline residue can be readily accommo- 
dated in the structure. It is interesting that both polyglycine and polypro- 
line are found to take up a configuration similar to the collagen chain (Rich 
and Crick, 1955; Cowan and McGavin, 1955). Thus a close similarity is 
to be expected between the structures of collagen and elastin, although so 
far there is no evidence that every third residue in elastin is glycine or 
regarding the relative disposition of glycine and proline residues. However, 
it is possible that the number of residues per turn (which is very nearly 34 
per turn in the minor helix) may not be exactly the same in elastin. 


4. RELATIONSHIP BETWEEN ELASTIN AND CHEMICALLY TREATED 
COLLAGEN 


Recently the authors have made a study of the effect of various inorganic 
salt solutions like those of nickel nitrate and calcium chloride on collagen. 
Collagen contracts by treatment with these reagents even at room tempera- 
ture. The progressive transformation of the structure brought about by 
these reagents has been followed by recording the X-ray patterns after 
treatment for various intervals of time (unpublished). In the case of both 
reagents, it is observed that there occurs first a disorientation of the struc- 
ture and then both the 2:86 A arc and 11 Ato 12 A equatorial spots simul- 
taneously disappear. Just after this stage, the pattern of the treated fibres 
shows a remarkable similarity to that of elastin. Patterns of the calcium 
chloride treated collagen and of normal elastin recorded on a flat film at 
3cm. are given in Figs.5 and 6. It will be noticed that the patterns are 
closely similar and that one could be easily mistaken for the other. The 
occurrence of a ring at about 9 A to 10A at the outer edge of the central 
halo in the pattern of the treated collagen specimen (exactly as with elastin) 
is particularly noteworthy. 


In the nickel nitrate treatment also, the same similarity is observable 
although there is a general background intensity probably due to the scat- 
tering and fluorescent X-rays from nickel in the pattern of this specimen. 
The ring at 9A to 10A has not been observed in the nickel nitrate treated 
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specimens if the concentration of the solution is large (100 per cent w/v), 
but the pattern given by them is closely similar to that of thermally shrunk 
collagen. If a dilute solution of nickel nitrate (30 per cent w/v) is used, 
and the collagen fibre is treated for a long time (several days), then the 
X-ray pattern of the fibre is very similar to that of the calcium chloride 
treated specimen, with a ring at about 9 A at the outer edge of the 
central halo. 


All these specimens also give the 2-2 A diffuse ring if the patterns are 
recorded in a 3cm. cylindrical camera. Usually the fibres begin to contract 
at the stage when the 2°86 Aand 11 A to 12A spots disappear, and they 
exhibit rubber-like elasticity at this stage. 


Further treatment with concentrated nickel nitrate solution causes the 
disappearance of the halo sorrounding the central spot and the pattern shows 
only a diffuse ring at 4-4 A with very little intensity inside. This is identical with 
the pattern of drastically shrunk rat tail tendon recorded by Astbury (1944) 
and of gelatin sol (high temperature form of gelatin) recorded by Katz et al. 
(1931). Elastin, when treated with nickel nitrate solution for a long time, 
also gives a pattern identical with the above pattern (see Figs. 7 and 8). It 
may be mentioned that the collagen fibre gives the usual pattern, if it is washed 
in distilled water, whatever be the stage to which it had been modified by 
the chemical treatment. 


These studies with chemically treated collagen confirm the ideas put 
forward in the previous sections. The fact that the treated collagen fibres 
can be readily restored to their normal condition by removing the salt means 
that the modification induced by the chemical treatment must be very mild. 
Since the X-ray pattern and physical behaviour of the specimen, treated 
under appropriate conditions, are very close to those of elastin, it is clear that 
elastin itself must be very similar to collagen. Further, it is observed that not 
only can collagen be modified further beyond this stage, when its X-ray 
pattern is similar to that of elastin, to that corresponding to hot gelatin, but 
elastin also exhibits such a change when treated with nickel nitrate. In 
fact, drastic thermal treatment also produces such a change in elastin, exactly 
as with drastically thermally shrunk collagen. If elastin is autoclaved and 
immediately the X-ray pattern is recorded, the pattern is identical with that 
of the drastically shrunk collagen, both showing only a fairly sharp ring at 
4-5A, with practically clear region within. This photograph is not repro- 
duced, but is very similar to Fig. 7(a@) of Astbury’s paper (1940) of rat tail 
tendon contracted to about 17 percent. of its initial length. 
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Similarly, it has been observed that if elastin is kept wet and the X-ray 
pattern taken, the pattern is very similar to that of wet gelatin. Thus, not 
only is it possible to produce a material from collagen which is very similar 
to elastin, but the further changes produced in both are also very similar. 


All these clearly show that elastin in fact belongs to the collagen family of 
proteins. 


5. OTHER OBSERVATIONS ON COLLAGEN AND ELASTIN 


When the above studies were nearing completion, an interesting article 
by Burton ef al. (1955) appeared, in which an apparent transformation of 
collagen fibres into a structureless elastin-like material was observed, when 
collagen was treated with some chemical reagents like borate buffer solu- 
tions. However, these authors have not observed the re-formation of the 
collagen, when the chemical is washed off. 


The similarity of shrunk collagen and elastin has also been observed 
by Banga (1953) who found that collagen is attacked and solubilised by 
elastase if it has been shrunk thermally or by treatment with phosphotung- 
stic acid, but is unattacked in the normal state. 


In view of the evidences from different sources, it appears to be very 
likely that proteins with properties and composition intermediate between 
those of collagen and elastin could also occur in biological systems (cf. 
Burton et al., 1955). It would be interesting to examine if this is so. The 
X-ray patterns of immature and developing collagenous tissue recorded by 
North ef al. (1953) appear to support this idea. With increase of age of rats, 
the pattern of the tail tendon is found to change from the elastin-shrunk 
collagen type to the normal collagen type. So also the infra-red data of 
developing tendon given by Jackson et al. (1953) show that the NH stretch- 
ing frequency (3310 cm.~*) of the developing tendon is significantly lower than 
that of adult tendon (3330 cm.—') and agrees with the value observed for hot 
gelatin (Robinson, 1953). It would be interesting to examine if elastin 
itself shows the lower frequency since the present studies indicate the similar- 
ity of the structure of elastin to degraded and less developed forms of collagen. 
So far no infra-red data on elastin appear to have been reported. 


Thus, elastin is probably a precursor of collagen in biological systems 
and collagen may be formed by the incorporation of polar residues and the 
resulting formation of cross-linkages. Burton et al. (1955) have on the other 
hand suggested that elastin is produced in vivo by degradation of collagen, 
although some of the co-authors (Hall et a/., 1955) themselves are of the view 
that this is probably less reasonable biologically than the reverse process. 
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6. SUMMARY 


The paper deals with the structure of elastin and its relationship to 
collagen. It is shown that, on the molecular level, elastin is analogous 
to collagen and is built up on a triple-chain helical arrangement of amino- 
acid residues. It is possible to obtain a material yielding an X-ray diffrac- 
tion pattern very similar to elastin, by treating collagen with calcium chlo- 
ride solution. The relationship between elastin and collagen in biological 
systems is also discussed. 
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DESCRIPTION OF PLATES 


II 


All the photographs have been taken with CuKa radiation in a crylindrical camera (R = 3 cm.) 


[Il 


X-ray diffraction pattern of calcium chloride treated collagen (Flat film, d = 3 cm,). 
X-ray pattern of elastin (Flat film, d = 3 cm.). 
X-ray pattern of nickel nitrate treated collagen (Flat film, d = 5cm.). 


Fic. 1. Elastin. 

Fic. 2. Collagen fibre. 

Fic. 3. Thermally shrunk collagen. 
Fic. 4. Gelatin. 

Fic. 5. 

Fic. 6. 

Fic. 7. 

Fic. 8. 


X-ray pattern of nickel nitrate treated elastin (Flat film, d = Scm.). 
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A NOTE ON CASCADE THEORY WITH 
IONISATION LOSS 
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(Department of Physics, University of Madras) 


Received January 3, 1957 
ABSTRACT 


The cascade theory of cosmic ray showers including ionisation loss 
is dealt with on the basis of the new approach suggested by us in an earlier 
contribution and an explicit Mellin transform solution is obtained for 
the mean number of particles produced in an infinite thickness of matter. 


Some years ago, Bhabha and Chakrabarthy (1943) obtained and discussed 
the equations of the cascade theory expressing the mean behaviour of the 
particles when collision loss is also taken into account. If f(E; 1) dE and 
g(E; ‘)dE are the mean numbers in the energy interval dE of electrons 
and photons respectively and R(E’\E) and p(E’|E) are the cross-sections 
for radiation and pair-production respectively, the equations are 


+ de’ + p 
= — (E31) dE 
+ RE — BIE) dE’ 2) 
with the initial conditions 
f(E; 0) = 6(E — E,) (3) 
where 5 is the Dirac delta function and 
g(E;0)=0 (4) 


corresponding to a shower excited by a single electron of energy Ey. In 
conformity with the theory of ‘product densities’* (which is usually useful 
for obtaining higher moments of the distribution) we recognise f and g as 


* For a formulation of the theory of product densities, see Ramakrishnan (1950). 
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— f(E; 1) J R(E’|E) dE’ + R(E|E) dE 
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the product densities of degree one of electrons and photons at ¢. The 
collision loss in traversing unit thickness of matter is assumed to be a con- 
stant, 8. Defining the Mellin transforms of f(E; f) and g (E; 1) as 


=f 1) (5) 


q(si) =) BME 6) 


according to Bhabha and Chakrabarthy, equations (1) and (2) can be reduced 
to 


= — Agp(s;t) + Beq(s3) —(s— (7) 
2 = C, p(s; t) — Dq(s; t) (8) 


As, Bs, C, and D are given by (see Bhabha and Chakrabarthy, 1943) 
4 


J 


We note that (7) is a difference-differential equation where the difference 
relates to the complex variable s and hence cannot be solved by direct 
iteration. Moreover the dependence on ¢ makes the solution much more 
difficult. A simple question that can be asked is: What is the spectrum at 
t= oo? Unfortunately the product density functions f) and g(E; 
tend to zero as t—>co (for E> 0) and thus nothing interesting can be 
obtained. However on the basis of the new approach suggested by us in 
a previous paper (1956) we can ask for the mean number of particles pro- 
duced in the entire shower, each of the particles having an energy greater 
than E, at the point of its production. We shall show that explicit Mellin 
transform solution for the mean number of particles in this case can indeed 
be obtained using standard mathematical techniques. 


| at 
an 
val 
in 
lc 
el 
[ 
V 
D — a 
9 6 


) 


' 


A Note on Cascade Theory with Ionisation Loss 135 


Let F (E; t) dEdt and G(E; t) dEdt be the mean numbers of elcetrons 
and photons produced between f and ¢ + dt respectively with their energies 
at the points of their production lying between E and E+ dE. F(E; ¢) 
and G(E; f) are product densities of degree one with respect to the two 
variables E and t. Using elementary probability arguments as we have done 
in our previous paper, we obtain 


F(E; 1) = t) p (E\B’) dE’ (10) 
G(E; 1) f(E'; 1) R(E’ — BIB’) dE’ (11) 


Note that (10) and (11) are valid whether or not we take ionisation 
loss into account. In fact, the dependence of F and G on ionisation loss 
enters through f and g. 

(10) and (11) are amenable to the standard Mellin transform technique. 


Defining P(s; 1) and Q(s; 1) as the Mellin transforms with respect to E, 
we find 


P(s; = | F(E; 1) = By q(s; 1) (12) 
Q(s; 1) =f G(B; 1) EE = C, P(s; 1). (13) 


The mean numbers of electrons and photons produced in the shower between 
0 and ¢ with the ‘ primitive’ energies (energy at the point of production) 
lying above E, is given by 


{N (Eg; = fdr f dE P(s;7)E-ds (14) 
0 Ee go — ico 


That these expressions tend to finite limits as ¢ + oo is obvious from 
a physical point of view since the number of particles with non-zero 
primitive energy by an initial particle of finite energy must be finite. 
Inverting the order of integrations we write the limits as 


{N (Ee; 00)} = 


ico 
= § q(s; 7) dr (16) 
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Thus the problem of obtaining the mean numbers reduces to that of solving 
for «,(s) and «,(s) where 


(s) = p(s; dr 


ka(s) =fq(s; 7) dr (18) 


We next proceed to show that though we cannot obtain explicit 
solutions of p(s; t) and q(s; t) from (7) and (8), we can solve for «,(s) 
and x. (s) quite easily. Integrating both sides of equations (7) and (8) with 
respect to ¢ from 0 to co and observing that p(s; 0) = E,S-! and qg(s; 0) =0 
we obtain 


— = (s) + (s) — (s — 1) Bx, (s — 1) (19) 
0 = Cgx, (s) — Dr, (s) (20) 
Eliminating «.(s), we have 


= — Age, (9) + Ba (s— QD 


+ - Asa) » 


(22 
K (8) = 


This recurrence relation can be continued and we have 


m= Asi) Asia): (= Asir) 


r=0 
Bms(s+ (s+ m) 


From (22), we note «,(s + m) for very large m is given by 


+m) (24) 
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so that the last term in (23) is negligible if m is sufficiently large. Hence we 
can write formally 

“= 4 (BD)'s 


r= 


(25) 


We shall show that the infinite series on the R.H.S. of (25) does converge 
absolutely so that (25) is the solution for «, (s). 


Writing 
Ky (Ss) = m) (26) 
we find 
¢ (m) DF 2 
d(m+1) Ey CsimBsim — DAsim 
For very large values of m, Cs,m and Bg,m are negligible and hence we write 
o(m) B .s+m 
For large m 
4 d 
Aum (29) 


The asymptotic expression for d/dz log \z can be derived by standard tech- 
niques [see for example Jeffreys and Jeffreys (1950), p. 466] as 


d 
08 = logz + (30) 
Thus 
d(m+1)| ~ Ey (3 «) (s + m) 
so that 
d(m + 1) 


tends to infinity as m—» oo. Hence the series X ¢(m) converges absolutely. 
Substituting (25) into (15) and (16), we find 


{N (Ee; D—1 I) 
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— 


i=1 


co ico 
1 Eo s-1 Bs 
(Ec; = (BD)™ f a 
m=0 


— 
(33) 


We note that the series in (32) and (33) contains inverse powers of 
as contrasted with the assumed expansion in terms of positive powers of § 
obtained by Bhabha and Chakrabarthy (1943, 1948). It is interesting t6 
note that in the case of a simple multiplicative model discussed in connet 
tion with the theory of multiple points by us elsewhere, the solution fog 
the mean number is also a series of inverse powers of f. 


We are thankful to Dr. V. K. Balachandran for discussions regarding 
the validity of our iterative procedure. 


One of us (S. K. S.) is grateful to the University of Madras for @ 
Government of India Scholarship. 
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